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1  Introduction 


Waves  and  particles  are  the  main  constituents  of  the  world  but  waves  are  our  main  method 
of  communication  and  detection,  whether  biological  or  technological.  Historically,  the  study 
of  waves  has  been  separated  into  the  stationary  (standing  wave)  and  nonstationary  (pulse) 
case.  Pulses  are  of  fundamental  consideration  in  radar,  sonar,  acoustics,  fiber  optics,  among 
many  other  areas.  A  pulse  has  been  called  by  many  terms:  transient,  wave  group,  progressive 
wave,  wave  packet,  nonrecurrent  wave,  traveling  wave,  non-periodic  wave  have  all  been  used 
to  describe  basically  the  same  thing.  Because  one  of  the  main  properties  of  pulse  propagation 
is  dispersion,  that  is,  the  fact  that  different  frequencies  propagate  with  different  speeds, 
time-frequency  analysis  offers  an  ideal  approach  toward  their  study  [3,  4],  Over  the  past  two 
decades  or  so  there  has  been  substantial  progress  in  the  development  of  new  methods  for 
analyzing  a  signal  jointly  in  time  and  frequency  and  this  has  led  to  new  results  in  the  basic 
nature  of  pulses. 

In  this  report  we  deal  with  pulses  and  one  of  the  main  aims  is  to  give  a  simple  view 
of  pulse  propagation,  an  approach  that  we  call  “a  local  view”  which  makes  pulses  easier 
to  understand  and  deal  with  intuitively.  We  have  devised  a  new  method  to  study  pulse 
propagation  in  dispersive  media  and  using  this  approach  we  have  derived  exact  expressions 
for  many  physical  quantities  such  as  the  spreading  of  a  pulse,  the  conditions  as  to  when 
a  pulse  contracts,  the  contraction  time,  and  other  important  physical  quantities.  We  also 
consider  higher  order  dispersion,  which  is  when  there  are  terms  on  the  dispersion  relation 
that  are  higher  than  quadratic,  a  concept  that  will  be  explained  subsequently. 

We  list  the  main  ideas,  results,  and  methods  presented  in  this  report: 

1.  We  describe  a  method  to  calculate  exact  moments  of  a  propagating  pulse  without 
having  to  calculate  the  pulse  itself. 

2.  We  discuss  a  local  view  of  pulse  propagation  that  explains  in  simple  terms  many  of  the 
curious  results  of  pulse  propagation. 

3.  We  have  introduced  the  concept  of  covariance  for  pulses  which  gives  a  clear  picture  of 
the  physical  situation  and  clarifies  many  of  the  mathematical  results  obtained. 

4.  We  have  studied  higher  order  dispersion  using  the  results  of  the  method  of  point  1. 

5.  We  show  how  instantaneous  frequency  changes  from  place  to  place  for  a  propagating 
pulse. 
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6.  We  have  explicitly  given  results  for  two  standard  situations,  when  a  pulse  is  known  for 
all  space  at  a  given  time  and  when  a  pulse  is  generated  at  a  particular  place  for  all 
time.  Both  situations  can  be  understood  using  the  local  view  model. 

7.  We  give  many  exactly  solvable  examples. 

8.  We  show  that  the  Wigner  distribution  approach  to  pulse  propagation  is  very  fruitful 
both  from  a  physical  and  numerical  point  of  view. 

Organization  of  this  report:  We  have  attempted  to  make  this  report  basic  in  regards  to 
the  issues  presented  and  hence  we  discuss  some  of  the  fundamental  facts  regarding  pulse 
propagation,  wave  equations,  and  dispersion.  In  this  manner  the  report  may  be  read  by 
someone  who  is  not  particularly  familiar  with  pulse  propagation  and  dispersion. 

2  Wave  Equations  and  Dispersion 

The  concept  of  dispersion  comes  about  in  the  following  way[12,  17] 
space  wave  equation, 

1  d2u  d2u 
c 2  d t2  dx2 

and  examine  how  the  simplest  wave  propagates 

us(x,  t )  =  eikx~iu}t  (2) 

If  one  takes  an  instant  snapshot  of  Eq.  (2),  that  is  if  one  looks  at  it  for  a  fixed  time,  then 
what  one  has  is  a  spatial  wave  and  its  spatial  frequency  is  given  by  k.  On  the  other  hand  if 
we  stay  at  one  position  then  what  we  have  is  a  wave  in  time  and  u  is  the  frequency.  Therefore 
both  k  and  u  are  frequencies,  one  a  spatial  frequency  and  the  other  a  time  frequency.  In  Eq. 
(2)  one  can  chose  k  and  lv  independently.  However,  what  the  wave  equation  does  is  force  a 
relation  between  the  two  and  that  is  called  the  dispersion  relation.  It  is  also  important  to 
emphasize  that  by  putting  Eq.  (2)  into  a  wave  equation  we  are  not  trying  to  find  a  solution 
but  just  seeing  if  such  a  solution  is  possible  at  all.  Also,  we  point  out  that  the  importance 
of  considering  Eq.  (2)  is  not  only  that  it  is  physically  simple,  but  because  any  pulse  can  be 
decomposed  in  terms  of  it,  where  k  and  to  range  over  all  space.  Hence  knowing  what  it  does 
for  k  and  to  allows  one  to  find  the  general  solution.  This  was  one  of  the  great  achievements 
in  pulse  propagation  and  is  discussed  in  subsequent  sections. 

We  also  should  point  out  that  getting  the  wave  equations,  Eq.  (1),  was  one  of  the  major 
problems  of  the  18th  century.  It  was  obtained  by  D’Alembert  and  Euler  and  each  had  a 


.  Take  the  standard  free 

(1) 
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different  view  point  as  to  what  is  an  acceptable  solution.  It  started  a  major  debate  that 
lasted  50  years  and  led  to  many  new  ideas  including  Fourier  analysis. 

If  one  substitutes  Eq.  (2)  into  Eq.  (1)  one  gets  that  it  can  only  be  a  solution  if  there  is 
a  relation  between  k  and  u.  In  particular 

u2  =  c2k 2  (3) 

or 

uj  =  ±ck  (4) 

Now  we  see  that  one  cannot  chose  k  and  u  at  will,  we  must  have  that  relation  satisfied.  Let 
us  momentarily  choose  to  =  k  and  therefore 

us(x,t)  =  eik{x~c^  (5) 

and  we  see  that  for  this  case  the  wave  propagates  to  the  right.  Moreover  one  can  now  show 
that  if  we  have  any  pulse  at  t  —  0  given  by  u(x,  0),  then 


u(x,  t)  —  u(x  —  ct,  0)  (6) 

That  is.  it  propagates  without  distortion.  That  Eq.  (6)  follows  is  because  an  arbitrary  wave 
can  be  decomposed  into  simple  waves  and  each  travels  at  the  same  speed. 

We  also  point  out  that  one  can  have  a  different  perspective  and  write 


us(x,t)  =  eiujWc~V 


and  we  will  discuss  this  later. 

Now  consider  the  beam  equation 

1  d2u  d  Au 
a2  dt2  dxA 

Again  putting  in  the  simple  solution  one  obtains  that 

a j2  =  a2k 4 


and  again  let  use  choose  u>  =  ak2.  Now 


us(x,  t )  =  e 


ik—iak2t 


_  ^ik(x—akt) 


(7) 


(8) 

(9) 


(10) 

(11) 


One  still  has  the  possibility  of  a  simple  wave  but  now  the  velocity  is  ak ,  that  is  the  velocity 
depends  on  k  and  therefore  if  we  had  the  sum  of  two  waves  they  would  not  move  in  unison. 
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Hence  if  we  have  a  pulse  that  is  the  sum  of  simple  pulses  there  is  no  one  velocity  associated 
with  the  whole  pulse  and  moreover  it  is  not  the  case  that  one  can  find  a  velocity  so  that 
Eq.  (6)  is  true.  There  are  other  velocities  that  can  be  associated  with  the  pulse  and  that  is 
something  we  will  discuss  in  detail  later. 

Generally  speaking  suppose  we  have  an  arbitrary  wave  equation,  put  the  simple  solution 


in  and  obtain  a  relation 

uj  =  W(k ) 

(12) 

then 

us(x,t)  =elk(x 

(13) 

and  we  see  that  a  simple  wave 

with  wave  number  k  propagates  with  a  velocity 

W(k) 

V*=  k 

(14) 

and  this  is  called  the  phase  velocity.  If  we  only  consider  simple  waves  this  would  be  it. 
However  if  we  consider  combination  of  pure  waves,  that  is  a  pulse,  these  considerations  are 
not  sufficient.  In  pulse  propagation  there  is  an  equal  and  perhaps  more  important  velocity, 
the  group  velocity,  as  we  will  explain  shortly. 


2.1  General  Solution 


Linear  partial  differential  equations  whose  solutions  give  wave  like  behavior  come  in  many 
varieties,  the  above  two  being  prime  examples.  Fortunately,  the  solution  to  all  of  them  can 
be  written  in  a  simple  form.  All  linear  wave  equations  with  constant  coefficients  may  be 
written  in  this  form 


N 


£ 

n= 0 


dnu 

'~dF 


M 


71=0 


dnu 

dxn 


(15) 


One  attempts  to  solve  it  by  substituting  elkx  lu3t  into  Eq.  (15)  to  give 

N  M 


(16) 


J2an(~iou)n  -J2bn(ik)n 

71=0  71=0 

One  can  solve  for  uj  in  terms  of  k 


u  =  W(k ) 


0 


(17) 

(18) 


or  solve  for  k  in  terms  of  u 

k  =  K(oj) 


(19) 
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Which  one  to  choose  depends  on  the  formulation  of  the  problem  and  specifically  on  the 
initial  conditions.  The  two  basic  initial  physical  situations,  that  we  call  case  A  and  case  B 
are 


Initially  Given  :  u(x,  0)  (Case  A)  (20) 

:u(Q,t)  (Case  B)  (21) 

The  first  case  is  when  we  have  the  spatial  wave  at  a  given  time,  and  the  second  when  we 
have  the  wave  at  a  given  position  for  all  time.  An  example  of  the  first  is  if  we  pluck  a  string 

and  let  go  at  time  zero.  An  example  of  the  second  is  if  we  are  at  a  fixed  position  and  create 

a  pulse,  for  example,  a  radar,  sonar,  or  fiber  optic  pulse. 

We  first  discuss  Case  A.  The  general  solution  is  [12,  17] 

u(x,t)  =  ~^=  f  S(k,  0)  eikx~iW{k)t  dk  (22) 

v27 t  J 

where  S(k,  0)  is  the  initial  spatial  spectrum  and  is  calculated  from  the  initial  pulse 

S(k,  0)  =  — L  [  u(x,  0)  e~ikx  dx  (23) 

V2vr  J 

and  W(k)  is  the  dispersion  relation  as  discussed  above.  The  general  solution  given  by  Eq. 
(22)  can  be  found  in  textbooks  [12,  17]. 

If  one  defines  the  time  dependent  spectrum  by 

S(k,  t)  =  S(k,  0)  e~iW{-k)t  (24) 

then  u(x,t )  and  S(k,t)  form  Fourier  transform  pairs  between  x  and  k 

u(x,t )  =  —  f  S(k,t)  elkx  dk 

v2vr  J 

S(k,t)  =  _  f  u(x,  t)  e~lkx  dx 

v2vr  J 

It  is  crucial  to  note  that  u(x,t)  and  S(k,t)  form  Fourier  transform  pairs 
for  all  time.  Case  B  will  be  developed  in  the  chapter  called  Case  B. 

Important  Concepts 

There  are  a  number  of  important  concepts  that  arise  and  we  discuss  them  here  briefly. 
Modes.  As  we  have  seen  from  the  examples  above  there  can  be  many  solutions  to  the 
dispersion  relation  and  each  solution  is  called  a  mode.  For  example  for  the  beam  equation 
we  have  that 

W(k)  =  ak 2  W(k)  =  - ak 2  (27) 


(25) 

(26) 

between  k  and  x 
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and  each  solution  is  called  a  mode.  The  general  solution  is  then  the  sum  of  the  modes.  In 
this  report  we  will  be  discussing  one  mode  at  a  time. 

Reality  and  attenuation.  If  for  a  particular  mode  W ( k )  is  real,  then  there  is  no  attenuation. 
The  reason  for  that  is  that  suppose  W ( k )  was  complex 

W(k)  =  WR(k)  -  iW^k) 

then  we  would  have 

^ ikx—iW(k)t  ^ikx—i(WR(k)—iWi(k))t 

—  e~Wi(k)t  ^ikx-iWR(k)t 

and  we  see  that  because  of  the  term  e~w,('k)t  the  wave  would  dye  out.  Therefore,  depending 
on  whether  W{k )  is  complex  or  real  we  will  have  damping  or  not.  In  this  report  we  consider 
the  case  where  we  have  no  damping.  However  we  emphasize  that  the  damped  case  is  very 
important.  See  the  “Future  Work”  section  of  this  report. 

Group  velocity.  A  central  idea  in  the  study  of  pulse  propagation  is  the  group  velocity,  v(k), 
which  is  defined  by 

v(k)  =  W\k )  (31) 

There  are  many  plausible  arguments  that  have  been  given  in  the  literature  for  calling  this 
quantity  the  group  velocity  and  most  books  give  a  plausible  argument  for  defining  the  above 
as  a  “velocity” .  In  a  subsequent  section  we  will  give  new  relations  for  a  propagating  pulse 
that  we  think  gives  a  very  clear  picture  as  to  the  physical  meaning  of  v{k )  is  and  how  it  is 
related  to  the  propagation  of  the  center  of  mass  of  the  pulse.  The  word  group  is  somewhat  of 
a  misnomer  but  it  comes  about  from  the  original  derivation  as  historically  given  and  indicates 
how  a  narrowband  set  of  waves  (a  group  of  waves)  centered  about  the  spatial  frequency  given 
by  k  propagates. 

Transit  time.  In  the  chapter  Case  B  we  will  study  the  time  properties  of  a  pulse  at  a  fixed 
position.  We  will  see  that  the  natural  quantity  that  appears  there  is 

t(u)  =  K'(pj)  (32) 

We  will  see  that  it  is  related  to  the  amount  of  time  it  takes  for  a  frequency  to  travel  a  unit 
distance.  We  call  it  the  frequency  transit  time. 

Structural/ geometrical  and  Media  dispersion.  As  we  have  seen  above  for  the  free  space 
wave  equation  there  is  no  dispersion.  But  that  there  is  dispersion  for  the  beam  equation. 
However,  if  one  imposes  boundary  conditions  on  the  free  space  wave  equations,  such  as  in 
a  waveguide,  then  dispersion  will  occur  because  the  imposition  of  the  boundary  conditions 


(29) 

(30) 
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generally  forces  a  relationship  between  u>  and  k.  This  type  of  dispersion  is  called  geometric 
or  structural  dispersion  in  contrast  to  media  dispersion  which  comes  from  the  wave  equation, 
such  as  the  beam  equation.  The  reason  one  has  dispersion  in  such  a  case  is  because  there 
is  reinforcement,  and  cancellations  of  waves  depend  on  the  geometry  and  therefore  certain 
modes  die  out  and  others  survive. 


2.2  Asymptotic  Solution 


One  of  the  fundamental  tools  for  linear  wave  propagation  has  been  the  so  called  asymptotic 
solution,  first  obtained  by  Kelvin  using  the  method  of  stationary  phase.  We  give  it  here 
because  most  of  our  work  will  be  to  obtain  exact  equations,  but  it  will  be  interesting  to 
contrast  with  the  traditional  asymptotic  solution.  We  do  not  give  the  derivation  here  as  it 
can  be  found  in  most  books.  One  takes  the  dispersion  relation  for  each  mode  and  sets 

W'(k)  =  J  (33) 

and  solves  it  for  k  which  we  call  ks 

ks  =  ks(x/t)  (34) 

These  are  the  so-called  stationary  points.  The  asymptotic  solution  is  then  [12,  17] 


ua(x,t)  ~  S(ks,  0)* 


tW"(ks 


^iksx— iW{ks)t— i7rsgn  W"  /  4 


(35) 


u, 


(  M\  S(ks,  0)  ikax—iW(ka)t—msgnW" /A 


=  Aa{x,t)ei^t) 

For  later  convenience  we  define  the  amplitude  and  phase 

’  y/tW"(ka) 

<pa{x,  t)  =  ksx  —  W ( ks)t  —  insgn  W" / 4 

A  simple  derivation  and  application  of  this  to  filtered  signals  was  given  by  Cohen  [9]. 


(36) 

(37) 


(38) 

(39) 


2.3  Pulse  Propagation  As  a  Filtering  Operation 

The  characteristic  of  pulse  propagation  is  that  the  spatial  spectrum  has  the  form 

S(k,t)  =  S(k,  0)  e~lW{k)t  (40) 
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It  is  a  product  of  the  initial  spatial  spectrum  multiplied  by  But  this  is  standard 

filtering  and  of  course  in  the  time  domain  it  is  the  convolution  of  the  initial  pulse  with  the 
function  that  produces  the  spatial  spectrum  e~lVV  ^ .  If  we  define 

w(x,t)=  1  f  e~iW{k)teikxdk 

v2n  J 

then  u(x,t )  is  the  convolution  of  u(x,  0)  with  w(x,t).  In  particular, 

(41) 

u(x,t)  =  Ju(x-x',oMx',tW 

(42) 

=  f  u(x',  0)w(x  —  x',  t)dx' 

(43) 

and  of  course  this  is  the  classical  expression  of  the  Green’s  function  approach.  The  reason  we 
mention  it  here  is  that  one  can  write  a  fast  and  simple  program  to  calculate  u(x,t)  because 
convolution  can  be  done  by  the  FFT  algorithm.  Dr.  Patrick  Loughlin  [13]  has  written  a 
simple  program  using  this  ideas  that  we  have  found  very  useful  to  study  the  issues  discussed 
in  the  report. 


3  Exact  Moments  of  a  Pulse 


We  define  the  moments  of  the  pulse  in  the  standard  way 

( xn)t=  j  xn  \u(x,  t)\2  dx  (44) 

and  our  aim  is  to  show  that  one  can  calculate  the  exact  moments  even  though  we  do  not 
solve  for  u(x,t).  There  are  three  reasons  for  this: 

•  First,  it  is  always  the  case  that  one  can  calculate  these  moments  from  the  spectrum, 
that  is 


( xn)t  =  J  x11  \u(x,  t)\2  dx 

=  [  S*(k,t)XnS(k,t)dk 


where  X  is  the  position  operator  in  the  k  representation 

-4 


(45) 

(46) 


(47) 


•  Secondly,  and  what  potentially  offers  a  simple  method  of  calculating  moments,  is  the 
fact  that 

S(k,  t)  =  S(k,  0)  e~iW(k)t  (48) 

and  therefore  the  differentiation  indicated  by  Xn  S(k,t)  can  actually  be  carried  on. 


•  Third,  once  the  differentiation  is  done,  the  resulting  quantities  are  expressible  in  terms 
of  polynomials  in  time  and  S(k,  0)  which  in  turn  means  we  can  express  things  in  terms 
of  the  initial  quantities. 

We  now  give  the  results  for  some  important  low  order  moments  and  discuss  the  physical 
meaning.  The  derivations  are  given  in  the  appendix. 

3.1  The  Mean 

The  first  conditional  moment  is 

( x)t=  [  x\u(x,t)\2dx 

S*(k,t)XS(k,t)  dk 

and  evaluates  to 


(49) 

(50) 


(x)t  =  (x)o+  vt  (51) 

where 

V  =  J  v(k)  |S(fc,0)|2dfc 
d\V(k) 

v(k)  =  sr 

We  now  discuss  the  physical  meaning  of  Eq.  (51): 

a)  First,  we  point  and  emphasize  that  the  usual  definition  of  group  velocity,  Eq.  (53), 
appears  in  a  natural  way  as  can  be  seen  from  the  derivation  in  the  appendix.  It  does  not 
have  to  be  imposed  in  any  way. 

b)  The  center  of  mass  travels  with  a  constant  velocity,  V,  for  all  wave  equations  and  all 
situations,  and  all  starting  conditions. 

c)  The  velocity,  V,  depends  of  course  on  the  dispersion  relation.  However  it  also  depends 
on  the  initial  spectrum.  Therefore  how  one  starts  the  pulse  effects  the  propagation  of  the 
center  of  mass! 

d)  Also,  V  depends  only  on  the  magnitude  of  the  initial  spectrum. 

e)  Note  that  V  can  be  considered  an  average  with  the  weighing  function  | S(k,  0)|2.  It  is 
the  average  of  the  group  velocity.  Therefore  v(k)  can  be  thought  of  as  the  velocity  for  each 
value  of  k  and  we  will  sometimes  write 


{v) 


J  v{k)  \S(k,  0)|2  dk  =  V 


(54) 
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3.2  Second  Moment 


and  evaluates  to 


(x2)t=  I  x  2\u(x,  t)\2  dx 

=  I  S*(k,t)X2  S(k,t)dk 

( x2  )t  =  ( x2  )o  +  t  ( vX  +  Xv  )o  +  t2  ( V 2  ) 


Where 


[vX  +  Xv)t  —  /  S*(k,t )  +  Aft7(A:)]  S(k,t)dk 


=  I  S*(k,t ) 
=  i  [  S*(k,t) 


n  s .  d  .  d 

v{k)'m+'dkv{k) 


S(k,t )  dk 


2 v(k)  — +  v'(k)S(k,  t ) 


dk 


A  standard  notation  for  such  quantities  is  the  anticommutator 


v,  x}+  =  v{k)X  +  Xv(k) 


(55) 

(56) 

(57) 

(58) 

(59) 

(60) 

(61) 


Also  this  quantity  can  be  calculated  from  the  pulse  directly.  If  we  take 

K  =  (62) 

then  also 

(uA  +  A’u)t  =  J  u*(x,t )  [u(/C)x  +  xu(/C)]  t)  dx  (63) 

This  quantity  and  its  physical  interpretation  will  be  quite  important  and  we  emphasize  that 
it  comes  up  naturally  in  the  derivation.  Alternative  expressions  will  be  given  in  the  next 
section,  where  we  will  also  discuss  its  physical  meaning.  Also,  in  Eq.  (57),  ( v 2)  is  obtained 
from 

(v2)  =  I  v2(k)\S(k,0)\2  dk  (64) 


3.3  Spread 

The  spread  at  a  particular  time,  defined  in  the  usual  way  by 

=  f  S*(k,t)(X  -  {x)t)2  S(k,t)dk 
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(65) 

(66) 


works  out  to  be  remarkably  simple, 

°l\t  =  <7x|o  +  2  *  Coyxv\o  +  t2a2v  (67) 

where 

°l  =  (v2)-(v)2  (68) 

=  j{v{k)  -V)2\S(k,0)\2dk  (69) 

Cov^io  =  \{vX  +  Xv)Q  -  (v)o{x)0  (70) 

We  now  make  some  remarks  regarding  the  above  expression  for  the  spread: 

1.  Notice  that  o2  is  the  standard  deviation  of  the  group  velocity.  Why  should  that  come 
in?  Notice  that  it  is  calculated  at  the  initial  time. 

2.  As  is  the  quantity  ( vX  +  Xv  )o  basic  to  our  considerations  so  is  the  quantity  CovOT,|o; 
we  call  it  the  “covariance”  between  position  and  group  velocity  because  it  acts  very  much 
like  the  standard  covariance.  It  is  also  calculated  at  the  initial  time. 


3.4  Covariance 


The  covariance  plays  a  basic  role  in  pulse  propagation  and  we  now  discuss  it  in  more  detail. 
In  the  above  we  defined  it  at  time  zero  but  now  we  define  these  quantities  for  an  arbitrary 
time.  In  particular  we  define 


(xv)t  =  \(vX  +  Xv)t  (71) 

=  i  J  S*  (k,  t)  {vX  +  Xv)  S(k,  t )  dk  (72) 

If  the  spectrum  is  written  in  terms  of  amplitude  and  phase 

S(k,t)  =\S(k,t)\ei^k't)  (73) 

then,  substituting  Eq.(73)  into  Eq.  (71)  one  obtains  that 

(xv)t  =  -  J v(k)  — ^  \S(k,t)\2dk  (74) 

Consider  how  (xv)t  changes  in  time  for  a  pulse.  The  time  dependent  spatial  spectrum  is 
given  by 


S(k,t) 


|S(M)|e#(M) 

(75) 

S(k,  0)  e~iW{k)t 

(76) 

|  S(k,  0)1  j'KW)e-iW(kyt 

(77) 

|  s(k,  0)1  eimo)-iw(k)t 

(78) 
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and  hence  the  phase  evolves  as 


ip(k,  t )  =  -0( k ,  0)  —  W ( k)t 


Differentiating, 


or 


di/j(k,t)  d'ip(k,  0) 

dk  dk 

d^(k,t)  d't/j(k,0) 


dk 

Substituting  this  into  Eq.  (74)  we  have 


dk 


-  W\k)t 
v(k)t 


(xv)t  =  -  [  v(k)  — \S(k,t)\~dk 


=  ~  J  v(k) 

=  -  [  v(k) 


dk 

di/j(k,  0) 

dk 

di/j(k ,  0) 


—  v(k)t 


\S(k,  t)\ 2  dk 


and  hence 


Now  consider  the  covariance 


dk 


\S(k,t)\2  dk  +  /  v(k)v(k)t\S(k,t)\2  dk 


(xv)t  =  (xv)o  +  ( v2)t 


Covx„| t  =  {xv)t  -  (x)t(v)t 


First  note  the  following 


( v)t  =  j  v(k)\S(k,t)\2 dk 

=  J  v(k)\S(k,  0)|2  dk 

=  (v)o 

Therefore,  the  average  group  velocity  does  not  change  and  we  have, 


(79) 

(80) 
(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 

(89) 


CovOT|,  =  {xv)t  -  (x)t(v),  (90) 

=  (xv)0  +  {v  2)«-((i)0+  Vt)(v)  o  (91) 

=  (xv)o-  { x)o(v)0  +  *((  w2>  -  V(v)o)  (92) 

=  (xv)o-  (x)o{v)o  +t({ u2)  -  (v)l)  (93) 

Hence 

Covxv\t  =  CovCT|o  +  (J2vt  (94) 


12 


Since  the  coefficient  of  t  is  manifestly  positive  we  see  that  no  matter  how  negative  the  initial 
covariance  is,  the  covariance  must  eventually  turn  positive.  Why  should  that  be  so?  What 
are  the  implications  of  that?  We  will  see  that  with  the  model  we  subsequently  present  these 
questions  will  become  easy  to  answer. 

Now  consider  the  correlation  coefficient  between  x  and  v, 


Covxv\t  _ _ Covct|q  +  ta2 _ 

ax\tav\t  av  ^ j 0-2  0  +  2 1  Cov^io  +  t2a2 


(95) 


We  now  examine  how  it  behaves  for  large  and  small  times.  For  large  times  we  have 


P  xv\ t 


(Cov 

XV  |o  +tcrl 


t<7 2 


1  + 

1  + 

1  + 

1  + 

1  + 


Cov 


xv\0 


ta2v 

Cov^jo 
ter 2 


Cov 


xv  1 0 


tol 

Cov^jo 

ta2 


1  + 


^x\0  "h  2  t  Cov-rqo 
t2a2 


-1/2 


1  + 


+  2 1  Cov^io 


t2<72 


-1/2 


1  a2j0  +  2 1  Cov^io  3  (  cr2|0  +  2 1  Covxq0 


t2a2 


+ 


|o 


t2a2 


Cov^io  1  0*|o  _ 


2  3/4  Cov2 


xv  1 0 


ta2  2  t2a2  8  \  t2at 

2  -  rr2  -2 

'J^v.r7,ln  a 

+ 


Cov^io  3  Cov^^jQ  ^x\o^v 


ta2 


2  t2at 


2  CovL|o  -  ol\o<rl 
2  t2oi 


As  time  goes  to  infinity  we  have  that 


(96) 

(97) 

(98) 

(99) 
(100) 
(101) 


Pxv\t 


as 


oo 


(102) 


which  shows  that  at  infinity  there  is  perfect  correlation  between  group  velocity  and  position! 
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For  small  times 


P  xv\ t 


Cov^iq  +ta. 


av\Jal\o  +  2t  Covx„|o  +  t2tf 


Cov 


xv\0 


tat 


&v&x\o  \J ^  (2 1  Cov xv  \o  +  t2tf)/crl\0 


Cov 


xv\0 


to, 


2  r 


&v&x\0 


1  -  ]^(2tCavxv\Q  +  t2a2v)/al\Q  +  jj(2 1  Cov^o  +  t2al)'2/a. 


4 

i|0 


_  CovHo  ,  yv  -  Cov^|0 

i  L 

&v&x\0  &v&x\0 

°v  -  CoVL|0 

—  p  ln  +  t - ^ — 

r  xv\0  2 

<rv<rx\o<7x\0 

3.5  Expansion  and  Contraction  Time 

We  repeat  here  the  equation  for  spread 

<rl[t  =  al\o  +  2  <  Cov„|0  +  (V2 


(103) 

(104) 

(105) 

(106) 

(107) 

(108) 


(109) 


and  emphasize  that 


•  The  spread  has  only  a  linear  and  quadratic  term  in  time. 

•  The  coefficient  of  the  quadratic  term  is  manifestly  positive. 

•  This  is  not  an  approximation  but  is  exact. 


The  quadratic  term  will  eventually  dominate  since  its  coefficient  is  manifestly  positive 
and  hence  the  spread  will  always  become  infinite  at  infinite  time.  But  the  coefficient  of 
the  linear  term  may  be  negative  and  therefore  the  pulse  may  contract  if  the  linear  term 
dominates  the  quadratic  term  for  a  period  of  time.  We  now  obtain  that  period  of  time.  If 
we  want  smaller  width  than  the  original  width,  that  is 

°l\t  <  ^|o  (no) 

then  from  Eq.  (109)  we  have 

<t*|o  +  2 1  CovHo  +  t2(Jv  <  al\o  (m) 
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or 

2 1  Covxt)|0  +  t2cr2  <  0 

giving  that  it  will  happen  for  contraction  times,  tc ,  given  by 

Cov  'fit  i  n 

0  <  tc  <  -2 - 

Therefore  for  there  to  be  contraction  we  must  have 

Covx„|o . -  0 


(112) 


(113) 


(114) 


and  why  that  should  be  the  case  will  be  seen  simply  when  we  explain  our  model.  To  obtain 
the  time  when  the  pulse  achieves  a  minimum  width,  differentiate  cbL  with  respect  to  time 
and  set  equal  to  zero, 


d°l\t 

dt 


=  2  Cov^o  +  2  tal  =  0 


yielding 


tM  — 


C'Ov.j.qo  1 
= 


(115) 


(116) 


ai  2 

therefore  the  time  to  achieve  the  minimum  contraction  is  equal  to  the  time  to  go  from  the 
minimum  contraction  back  to  the  original  width.  To  determine  how  thin  a  pulse  can  get  we 
substitute  t  =  tM  in  Eq.  (109)  to  obtain 


® x\tM  Tr|0  (l  Pxv |o) 


(117) 


where  pxv\ 0  is  the  correlation  coefficient, 

2 

Pxv  |0 


CovLio 

^^10 


(118) 


We  now  give  a  physical  picture  explaining  the  above  results.  Consider  the  initial  pulse 
to  be  a  swarm  of  particles  where  each  particle  has  a  constant  velocity  and  where  the  swarm 
has  an  average  velocity  V.  We  take  V  to  be  positive,  that  is  the  center  of  mass  is  moving  to 
the  right.  Now  suppose  that  initially  the  particles  on  the  right  have  high  velocity  and  the 
particles  on  the  left  have  relatively  slower  velocity.  Since  everybody  travels  with  a  constant 
velocity  the  spread  will  be  higher  once  the  particles  start  to  move.  On  the  other  hand 
suppose  that  fast  particles  are  to  the  left  and  slow  particles  are  to  the  right.  Once  the  swarm 
starts  to  move  the  fast  particles  are  catching  up  and  the  distance  between  the  slow  and  fast 
particles  is  decreasing.  That  is  the  swarm  contracts.  Eventually  the  fast  particles  catch 
up  and  that  is  the  time  when  the  pulse  width  is  minimum.  Once  they  catch  up,  they  pass 
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them  and  the  distance  then  continues  to  increase.  The  reason  the  covariance  comes  in  such 
a  fundamental  way  is  that  it  tells  us  whether  the  fast  particles  are  to  the  right  or  left  of  the 
pulse.  If  the  covariance  is  positive  the  fast  particles  are  to  the  right  of  the  slow  particles  and 
if  it  is  negative  it  means  they  are  arranged  the  opposite  way. 

We  now  explain  why  the  contraction  time  also  depends  on  the  spread  of  initial  velocities. 
If  the  spread  is  very  large  then  that  means  there  is  a  big  difference  between  the  fast  and 
slow  particles  and  hence  the  spread  in  position  will  decrease  faster  since  the  fast  particles 
can  catch  up  quickly  to  the  slow  particles.  Conversely  if  the  velocity  spread  is  small  then  it 
will  take  much  longer  for  the  fast  particles  to  catch  up  to  the  slow  particles. 

Spread  of  Uncertainty  Product.  Using  the  fact  that  for  pulse  propagation 

°l\t  =  al\o  (119) 

the  uncertainty  product  evolves  as 

o%\t<Tk\t  =  Wilt  +  ZtCoVn,  +  t2al\al  |0  (120) 

=  +  [2£Covot  +  t2a2]a2k]0  (121) 

The  uncertainty  product  eventually  goes  to  infinity  but  can  also  momentarily  decrease. 
When  the  uncertainty  product  is  large  the  signal  is  said  to  be  asymptotic.  Thus  in  pulse 
propagation  all  signals  become  asymptotic. 

4  Local  Model  of  Pulse  Propagation 

We  now  give  the  details  of  a  simple  physical  model  that  gives  the  identical  equations  as  pulse 
propagation  and  from  which  all  the  quantities  can  be  understood  simply.  More  importantly 
we  argue  that  the  model  can  be  used  to  visualize  and  understand  pulse  propagation  in 
general.  Motion  of  objects  is  simpler  to  visualize  and  understand  than  waves  and  certainly 
much  of  the  vocabulary  of  waves  is  derived  from  the  desire  to  make  believe  waves  act  like 
particles.  We  present  a  simple  model  of  particle  motion  that  is  very  close  to  wave  motion 
and  which  hopefully  will  help  to  understand  wave  motion  phenomena.  The  reason  one  may 
want  to  develop  a  particle  view  of  wave  motion  is  manifold.  First,  it  may  give  us  a  clear 
and  more  intuitive  picture  of  what  is  going  on  with  wave  phenomena.  Second,  it  could 
lead  to  practical  numerical  schemes  because  instead  of  solving  the  wave  equation  we  can 
evolve  particles.  Third,  it  may  help  to  understand  non  linear  waves  which  are  particularly 
difficult  to  visualize.  In  addition,  as  we  will  show,  in  trying  to  understand  wave  propagation 
one  comes  across  certain  physical  quantities  that  appear,  mathematically,  like  “correlation” , 
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although  their  mathematical  properties  do  not  fully  match  those  of  classical  correlations.  By 
making  a  particle  model  we  may  be  able  to  understand  better  correlations  between  quantities 
as  they  appear  in  the  study  of  wave  propagation. 

Consider  the  simplest  motion  of  a  particle,  the  motion  with  a  constant  velocity, 

xt  —  xq  +  vt  (122) 

where  xt  is  the  position  at  time  t,  and  the  constant  velocity  is  V.  Now  suppose  we  have  a 
group  of  particles  and  the  group  is  characterized  by  a  density  at  the  initial  time.  If  all  the 
particles  have  the  same  velocity  it  is  clear  that  the  group  will  move  but  keep  its  shape.  Now 
suppose  the  velocity  of  each  particle  is  not  the  same  and  suppose  we  pick  the  velocity  of 
each  particle  depending  on  the  initial  position,  x0.  Hence  we  write 

v  =  v(x0)  (123) 

and  now  we  have 

x(t)  =  xq  +  v(x0)  t  (124) 

For  the  density  of  particles  we  will  use  p(x,t),  that  is 

p(x,  t )  =  the  density  of  particles  at  time  t  (125) 

Also  for  convenience  we  use  the  following  notation  for  the  density  at  time  zero, 

p(xo)  =  p(x,0)  (126) 

It  will  be  helpful  to  think  of  this  common  situation.  Suppose  we  have  a  large  number  of 
runners  at  the  starting  line  of  a  marathon  and  let  us  assume  they  take  up  a  city  block. 
We  could  place  the  runners  in  many  different  initial  configurations.  The  most  common 
arrangement  is  where  one  places  the  fast  runners  at  the  starting  line,  but  we  could  place 
them  at  the  end  of  the  block,  or  in  the  middle,  or  any  way  we  want.  Each  of  these  different 
ways  is  characterized  by  the  function  v(x0).  How  will  the  group  of  runners  behave  once  the 
gun  goes  off?  Clearly  if  the  fast  runners  are  at  the  starting  line  the  swarm  will  immediately 
spread.  But  if  the  fast  runners  are  at  the  back  at  the  start  of  the  race,  then  the  swarm  will 
first  contract  and  then  expand. 

Now  let  us  study  how  the  group  of  runners  behaves  after  the  starting  gun.  As  usual  we 
define  moments  as 

( xn)t=  j  xnp(x,t)dx  (127) 

However  we  do  not  have  to  know  p(x,  t )  explicitly  because  we  can  calculate  everything  from 
the  initial  density 

{xn)t=  j (x0  +  v(x0)  t)np(x0)dx  (128) 
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First  consider  the  center  of  mass.  One  takes  expectation  values  of  both  sides  of  Eq.  (124) 


(xt)  =  (x  o)  +  (v(x0))t 


(129) 


or 

(x)t={x)0  +  Vt  (130) 

where 

V  =  (v0)  =  J  v(x0)p(x0)dx0  (131) 

Eq.  (131)  shows  that  the  center  of  mass  moves  with  a  constant  velocity  which  is  the  average 
velocity  of  the  group  at  the  initial  time.  Now  consider  the  second  moment.  Square  Eq.  (124) 

xt  =  (x0  +  v0 1)2  =  xl  +  2  x0v0t  +  vlt2  (132) 

and  take  averages  of  both  sides,  giving 

( x2  )t  =  ( x2  )o  +  2(  xv  )0t  +  ( v 2  )ot2  (133) 


where 

{x0v0)  =  j  xov(xo)P(xo,0)dx0 
The  spread  of  the  swarm  is 

°l\t  =  (x)2t-(xt)2 

=  (x2)o  +  2(xv)o  t  +  (v2)0 1 2  -  ((x)o  +  {v)o tf 
=  ( X2)o  -  {x)l  +  2((xv)0  -  (x)o(v)o)t  +  (( v2)o  -  (■ v)l)t 2 


(134) 


(135) 

(136) 

(137) 


or, 

°l\t  =  ax\o  +  2  Covra|0t  +  v2]0t2  (138) 

where  cr^0  is  the  standard  deviation  of  the  initial  velocity 

^2|0  =  <^2>0  -  (v)l  (139) 

and  where  the  covariance  is 

Covra|0  =  (xv(x))o  -  (v)o(x)o  (140) 

It  is  the  covariance  between  the  initial  position  and  initial  velocity.  One  can  write  the  spread 
in  terms  of  the  correlation  coefficient, 

ax\t  =  ax\o  +  2pvx\o^x\o^v\ot  +  crv\0t  (141) 
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where  the  correlation  coefficient  is 


Pvx  |0 


^x\0  &v\0 


(142) 


With  one  exception,  the  coefficient  of  t2  in  Eq.  (141)  can  never  be  made  zero  and  therefore 
the  swarm  must  spread  to  infinity.  The  coefficient  of  t2  can  be  made  zero  only  if  we  take 
v0  =  c.  For  that  case  we  have  that  Covxvjt  and  cr^0  are  both  zero  and  hence  a2,t  =  cdj0  and 
the  swarm  will  not  spread.  While  it  is  certain  that  the  swarm  will  spread  eventually,  the 
swarm  can  contract  for  a  certain  time.  For  the  pulse  to  contract  we  must  have  Covxv  <  0 
which  gives 

0  <  t  <  _2  Co^bo  (143) 

av\0 

for  the  contraction  time. 

Let  us  also  calculate  the  covariance  as  a  function  of  time.  First,  we  have  that 


x(t)v(t)  =  (x„  +  vot)vo  =  Wo  +  v„ 4 


(144) 


Taking  averages  we  obtain 


The  covariance  at  time  t  is  then 


{xv)t  =  (xv)0  +  (vjj)t 


(145) 


Covvxjt  =  (xv)t  -  (x)t(v)t 


(146) 


which  gives 


CoVXV\t  C<OVa;a;|o  T  CT„|q  t 


We  see  that  the  covariance  must  increase.  For  the  correlation  coefficient  we  have 


Pxv\t 


Covx„|t  _  Cov^io  +  av{0t 

XxltXvlt  crv|o  \fo%o  +  2  Cov^iol  + 


(147) 


(148) 


The  covariance  and  correlation  coefficient  increase  because  as  the  swarm  spreads  the  fast 
particles  will  be  more  and  more  up  front  of  the  swarm  and  will  hence  be  highly  correlated 
with  the  relative  position  in  the  swarm.  Indeed  as  time  goes  to  infinity  the  correlation 
coefficient  goes  to  1, 

Pxv\t  — >  1  as  t  oo  (149) 

To  understand  these  equations  consider  the  case  where  the  fast  runners  are  lined  up  in 
the  front.  The  initial  covariance  is  positive,  which  indicates  that  fact.  At  the  start  the  fast 
runners  take  off,  the  slow  runners  at  the  back  do  not  keep  up  and  the  swarm  immediately 
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spreads.  However,  suppose  that  the  fast  runners  are  in  the  back:  the  covariance  is  then 
negative  which  means  that  the  slow  runners  are  np  front.  Now,  for  a  short  time  after  time  at 
the  start  of  the  race  the  fast  runners  are  catching  up  to  the  slow  runners  and  therefore  the 
swarm  contracts.  Also,  the  bigger  the  initial  difference  between  the  fast  and  slow  runners, 
the  longer  will  be  the  contraction  time.  In  addition  the  bigger  the  difference  in  speed  between 
the  fast  and  slow  runners  the  smaller  will  be  the  contraction  time,  because  the  fast  runners 
will  catch  up  more  quickly.  The  same  type  of  behavior  is  exhibited  by  pulse  propagation  in 
dispersive  media. 


5  Instantaneous  Frequency  and  Dispersion 

“Instantaneous  frequency”  is  one  of  the  most  fundamental  quantities;  it  is  frequency  as 
a  function  of  time.  While  it  is  intuitively  obvious  that  frequencies  do  change,  its  exact 
mathematical  expression  is  far  from  obvious.  Generally  speaking  instantaneous  frequency 
is  the  derivative  of  the  phase  and  the  question  has  been  how  do  we  determine  the  phase. 
This  is  an  old  and  important  problem  and  the  general  solution  is  usually  the  one  given  by 
Gabor.  Take  the  real  signal  and  delete  its  negative  frequencies  and  form  a  new  signal  just 
from  the  positive  ones.  That  results  in  a  complex  signal  and  the  derivative  of  the  phase  is 
then  the  instantaneous  frequency.  However  this  procedure  has  not  been  fully  investigated 
for  a  pulse  and  we  do  not  do  so  here  but  just  point  this  out.  Here  we  use  the  solution  to  the 
wave  equation  and  express  it  in  terms  of  amplitude  and  phase, 

u(x,t)  =  (150) 


to  define  instantaneous  frequency  as 


Vi(x,t)  =  ~—(p(x,t) 

and  spatial  local  frequency  by 

d 

H- M)  =  gTfiX.t) 


(151) 


(152) 


The  reason  for  the  differences  in  sign  is  because  of  the  basic  definition,  Eq.  (2). 

We  now  address  the  problem  as  to  how  does  the  instantaneous  frequency  vary  from  point 
to  point  as  a  pulse  propagates  and  how  does  dispersion  effect  that.  For  example  suppose 
a  sound  is  made  and  there  are  two  people  each  standing  at  50  and  100  feet  away.  Will 
they  hear  the  same  sound?  The  answer  is  yes  because  for  frequencies  that  humans  speak  at 
there  is  no  dispersion  in  air.  However  suppose  there  is  dispersion,  how  would  the  frequencies 
change?  That  is  how  does  the  instantaneous  frequency  of  a  pulse  depend  from  position  to 
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position?  We  have  not  been  able  to  do  this  exactly  but  we  have  done  it  approximately  using 
the  asymptotic  solution.  However,  for  many  examples  we  have  been  able  to  do  it  exactly  and 
hence  we  can  compare  our  approximate  solution  with  these  exact  solutions.  For  convenience 
we  repeat  here  the  asymptomatic  solution 


ua(x,t )  ~  S(ks,  0) 


tW"(ks 


gikgX—iW  (ks)t—ii rsgn  W”  /  4 


and  the  amplitude  and  phase  are 

A«M  =  \S(k„  0)lyL^-y 

tpa(x,  t)  =  ip(ks,  0)  +  ksx  -  W(ks)t  -  7Tsgn  W" / 4 


(153) 


(154) 

(155) 


where  the  value  of  ks  is  obtained  from  solving 


W'(ks)  =  x/t 


Differentiating  the  phase,  (pa(x,t),  we  have 

r) 

Ui(x,  t )  = 


dtp  dks 

dk. 

s  dt 

dip 

dks 

r  dtp 

-I-  x _ 1  _  t _ v  _ f 

dt  dk ,  dt 


+  W(ks 


dW(ks 
dk , 


w+w™ 


dk , 


w  +  w™ 


and  hence 

From  Eq.  (156)  we  have 

W"(ks)^f  =  -x/t2 

giving 

dks  x 

~dt  =  ~t2W"(ks ) 
W'(ks) 
tW"{ks ) 

w2(ks ) 

xW"{ks) 


(156) 

(157) 

(158) 

(159) 

(160) 

(161) 

(162) 

(163) 

(164) 

(165) 
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and  therefore  we  can  express  the  instantaneous  frequency  in  the  asymptotic  regime  in  a 
variety  of  ways, 


x  dtfi 
t2W"(ks)dks 
W'(ks)  d xfj 
tW"(ks)dks 
W'2{ks )  d. 0 
xW"(ks )  dks 


■  +  W(kg) 

(166) 

+  W(ks) 

(167) 

+  W{kg) 

(168) 

ft  should  not  be  concluded  from  these  expressions  that  the  instantaneous  frequency  varies 
in  any  specific  way.  For  example  by  looking  at  Eq.  (166)  one  may  think  that  is  a  linear 
function  of  x.  That  would  not  be  right  because  in  all  of  the  above  examples  ks  is  a  function 
of  x  and  t.  However  one  can  get  a  sense  of  that  if  we  expand  the  dispersion  relation  in  a 
power  series. 

Spatial  Local  Frequency.  We  now  consider  the  issue  of  spatial  instantaneous  frequency. 
We  define  it  by 

ki{x,t)  =  ^-ip{x,t)  (169) 

ox 

Differentiating  the  phase,  (pa(x,t),  we  have 


d 

ki(x,t)  =  —<pa(x,t) 
ox 


dili  dks  dks 

dg^+x^+k‘ 
_  dW{ks) 
dk.  dk 


dil> 
dks 
dill  dks 
dks  dx 


+  x 


X 

t— 

t 


S 

dks  , 

- ^  ^s 

dx 


tdW{ks)dks 

dks  dx 

dks  , 

- ^ 

dx 


+  ks 


But  from  Eq.  (156)  we  now  have  that 


W"{kg) 


dks 

dx 


1/t 


giving 

dks  1 
~dx  =  tW"(kg) 
W\ks) 
xW"{kg) 


(170) 

(171) 

(172) 

(173) 

(174) 

(175) 

(176) 

(177) 
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and  hence 


k(x,  t ) 


1  dip{k,  0) 
tW"(fcs)  dA;s  s 
W'%)  #(fc,  0) 
xW"(ks)  dks 


(178) 

(179) 


This  is  the  spatial  instantaneous  frequency. 

Translation  invariance.  We  point  out  that  while  the  exact  solution  is  translationly 
invariant  the  asymptotic  solution  is  not.  In  particular  suppose  we  define 


utr(x,  0)  =  u(x  —  Xq,  0)  (180) 

then  for  the  exact  solution 

Utr(x,  t)  =  u(x  —  Xq,  t)  (181) 

Therefore,  if  we  work  things  out  for  u(x,  t )  it  is  easy  to  take  all  formulas  and  change  them 
for  the  case  u{x  —  Xq ,t).  One  merely  substitutes  x  —  xq  for  x.  However  for  the  asymptotic 
solution  that  is  not  the  case,  but  we  can  obtain  a  formula  to  handle  the  situation.  In  general 
we  have  that 


Str(k,t)  =  eikx°S(k,  t) 

(182) 

and  therefore 

ua,tr(x,t)  =  elksX0ua(x,t) 

(183) 

The  phase  hence  changes  according  to 

Ta,tr(x,  t )  =  (pa(x,  t )  +  ksX0 

(184) 

and  the  instantaneous  frequency  is  therefore 

.  .  dks 

^ i,tr\X ,  t)  U>iyX,  t )  Xq 

(185) 

(  ,  xxo 

=  ^t)  +  t2w„{ks) 

(186) 

similarly 

dk 

ki,tr(x,t)  =  ki(x,t )  -  -Q^xo 

(187) 

=  k^t)-tw>'lk,) 

(188) 
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6  Higher  Order  Dispersion 


Historically  for  a  variety  of  reasons  dispersion  relation  which  are  at  most  quadratic  have 
been  intensively  studied.  Perhaps  the  main  reason  is  that  for  that  case  certain  integrals  can 
be  carried  out.  Higher  order  dispersion  is  when  there  are  terms  in  the  dispersion  relation 
that  are  higher  than  quadratic  and  it  is  important  to  understand  their  effect.  Now  in  our 
previous  development  we  have  shown  how  to  obtain  exact  relations  for  any  dispersion  relation 
and  therefore  we  should  be  able  to  get  exact  results  illustrating  higher  order  dispersion.  To 
obtain  a  sense  of  higher  order  dispersion  we  first  do  a  specific  example  and  then  see  if  with 
the  simple  example  there  are  qualitatively  different  results.  If  so  then  we  will  try  to  prove 
these  in  general.  The  example  we  use  is  where  we  add  a  cubic  term  to  the  dispersion  relation 
[10] 

W(k)  =  ck  +  q/c2/2  +  r]k3/3  (189) 

and  for  the  initial  pulse  we  take 


u(x,  0)  =  (a/V)1/4  e-ax2 /z+iP*2 /2+ikox 


(190) 


For  each  physical  quantity  of  interest  we  shall  give  the  result  and  discuss  the  effects  of  the 
cubic  dispersion  term. 

Velocity  of  the  center  of  mass.  One  obtains 


V 


c  +  jk  o  +  T) 
c  +  q/co  +  V 


k20  + 


kl  + 


2  a! 
a 2  +  (3 2 
2a 


(191) 

(192) 


Without  the  cubic  term,  it  is  seen  that  what  enters  are  c  and  k o,  and  that  is  reasonable 
because  k o  is  just  the  average  spatial  frequency  of  the  initial  pulse.  However,  adding  the 
cubic  term  introduces  fundamental  new  parameters,  namely  a  and  /3,  which  were  absent 
before. 

Spread  of  Group  Velocity.  The  spread  is 

1  ' 


<T„2  = 


7 

2a' 

1 

a1 


2  rnkp 


a 


rf 
+  V 


2k°  +  M 


T 


2 

a2  +  /3 2 


+  77^o  +  V2 


2  kn  + 


2a' 


a 


y  +  77  +  V2 


2k°  +  2^ 


(193) 

(194) 

(195) 


Covariance.  As  we  discussed  the  sign  of  the  covariance  is  crucial  in  determining  whether  a 
pulse  will  contract  and  for  how  long  the  contraction  will  go  on,  and  therefore  it  is  important 
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to  see  how  higher  order  dispersion  effects  that.  The  covariance  is  now  given  by 


P 


CovT1.  =  — 


a  L2 


7T  +  Vko 


(196) 


The  higher  order  results  now  depends  on  k0.  If  we  only  have  the  quadratic  term  then  the 
sign  of  the  covariance  is  only  determined  by  the  sign  of  P,  assuming  that  both  a  and  7  are 
positive.  But  now  we  have  another  possibility  because  the  sign  and  magnitude  of  ka  is  at 
our  disposal.  In  fact,  we  can  make  the  covariance  less  than  zero  even  if  we  take  all  other 
parameters  positive.  This  will  be  the  case  when 


Cov  xv  <  0 


for 


k0<  ~ 


a  7 


V  IP  2_ 

Spread.  Putting  the  above  results  into  Eq.  (67)  we  have  that  the  spread  is 


<rl\t  =  °l\n  +  2 1 


’x\0 


—  ax\o  +  2t 


—  -.2 
ax\0 


[- 

'7:  +  Vko) 

a 

V2  ) 

- 

^  +  Vko) 

a 

\2  J 

+  r 


7  ,  ivkp  rj_ 
2  a'  01'  a 


2  /  1 

2/Cn  + 


2  m,2 


+ 


a' 


7 


+  rnhi  +  v2  2/cq  + 


2a! 

1 


2a' 


1  +  2/3t(j  +  2 arjk0)  +  ( a 2  +  P2)t2  <  q2  +  2jrik02r]-  2 k£  + 


a2  +  P2 
2a 


Contraction.  Using  Eq.  (198)  the  pulse  will  contract  for 

0  <  t  <  T 


(197) 

(198) 

(199) 

(200) 

(201) 


where 


T  = 


/3a' 

a 


+  Vko ) 


72/2  +  777&0  +  rf  (2 k20  +  ^7) 


2  P 


q(a2  +  P2) 


1  + 


V2  i2ko  +  jg) 

(1  +  2  i]k0) 


-1 


(202) 

(203) 


If  one  compares  the  cubic  case  with  the  quadratic  case  we  see  that  an  important  new  signal 
parameter  comes  into  play,  namely  h0.  It  is  probably  the  case  that  this  parameter  would  also 
enter  into  the  quadratic  case  if  an  initial  pulse  was  taken  so  that  it  would  not  be  symmetrical 
about  k0.  Nonetheless  what  our  example  shows  is  that  there  is  a  wide  variety  of  possibilities 
to  control  the  contraction  and  rate  of  expansion  of  a  pulse.  It  is  important  to  appreciate  that 
the  contrasts  may  just  be  due  to  the  fact  that  we  have  used  a  Gaussian  example  and  that 
perhaps  for  the  quadratic  term  the  peculiarities  arise  just  because  of  that.  That  is,  perhaps 
the  reason  the  cubic  terms  effects  arise  is  artificial,  and  perhaps  they  really  are  there  in 
general  but  absent  from  this  examples.  Further  studies  are  needed  to  understand  this  fully 
but  what  the  above  shows  is  one  can  calculate  appropriate  quantities  exactly  for  the  sake  of 
comparison. 
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7  Wigner  Distribution  Approach 


We  now  discuss  how  the  Wigner  distribution  [18,  3,  4]  gives  considerable  insight  and  calcula- 
tional  advantage  into  pulse  propagation.  It  is  natural  that  it  should  be  the  case  since  the  aim 
of  time-frequency  analysis  is  to  show  how  frequencies  change  in  time.  The  spatial/spatial- 
frequency  Wigner  distribution  is  [11,  15] 

W(x,k,t)  — —  [  u*{x  —  hr,t)  u(x  +  hr,t)  e~lTk  dr  (204) 

27 r  J 

and  in  terms  of  the  spatial  spectrum  it  is 

W(x,  k,t)  —  —  [  S*(k  +  ±0,  t)  S(k  -  ±0,  t)  e~idx  dd  (205) 

27 T  J 

We  first  give  the  important  moments  in  a  general  way,  that  is  without  specifying  that  we 
are  dealing  with  pulse  propagation.  Before  giving  these  moments  we  repeat  here  the  basic 
definition  of  phases  and  amplitudes  for  both  the  signal  and  spectrum 


u(x,t )  =  A(x,t)el'p(x’^  (206) 

S(k,t )  =  _  f  u(x,t )  e~lkx  dx  (207) 

v27 t  J 

=  B(k1ty^t)  (208) 


The  Wigner  distribution  satisfies  the  marginals, 


J  W(x,  k,t)dk  =  \u(x,t)\2 

(209) 

J  W (. x ,  k,  t)dx  =  |  S(k,  t)  2 

(210) 

For  the  sake  of  clarity  in  this  section  we  will  not  normalize  the  moments  with  their  respective 
marginals. 

A  straightforward  calculation  yields 

( k)xt —  j  kW(x,k,t)dk 

(211) 

d  .  . 

=  ^,t) 

(212) 

which  is  the  instantaneous  spatial  frequency.  Also, 

(. x)kt  =  B(k,t )  J  xW(x,k,t)dx 

(213) 

=  -Av>(m) 

(214) 
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and  the  covariance  is 


(. kx)t  —  J  xkW(x,  k,t)dxdk  (215) 

-  [  xd(p^,t]\ u(x,t)\2dx  (216) 

J  OX 

=  _/4^M)|sftt)|2  (2l7) 

Now  let  us  specialize  to  the  case  of  pulse  propagation.  We  repeat  here  for  convenience 
the  basic  relevant  equations  that  characterize  pulse  propagation 


S(k,t)  =  S(k,  0)  e~iW{k)t 

(218) 

=  5(M)e#(M) 

(219) 

=  B(k,  o)e^(k^-iW(k)t 

(220) 

and 

B(k,  t )  =  B(k,  0) 

(221) 

tp(k,  t )  =  ip(k,  0)  —  W ( k)t 

(222) 

Therefore  using  Eq.  (222)  we  have 

d 

v x)k,t  =  -  <^(M) 

(223) 

=  -^(k,0)  +  W'(k)t 

(224) 

=  -^ip{k,0)  +  v(k)t 

(225) 

But  at  t  —  0 

IT 

j*- 

© 

J? 

o 

(226) 

and  therefore 

(x)k,t  =  (x)k,  o  +  v(k)t 

(227) 

This  is  an  interesting  result  because  it  predicts  exactly  the  picture  we 

developed  previously. 

Now  consider 

(kx)t  =  J  k  ^  \s(k,t)  |2 

(228) 

=  -/  k{^(k,o)-Ak)t)\s{k,o)\> 

(229) 

=  (kx)  o  +  (kv(k))0t 

(230) 
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which  agrees  with  our  previous  results.  The  important  point  to  observe  here  is  that  these 
expressions  and  results  came  out  naturally  using  the  Wigner  distribution. 

Evolution  of  the  Wigner  distribution  and  an  approximation.  Our  aim  is  to  express  the 
Wigner  distribution  at  time  t  in  terms  of  the  Wigner  distribution  at  time  t  =  0.  The 
reason  we  consider  this  problem  is  because  a  considerable  simplification  occurs.  For  pulse 
propagation  we  substitute  Eq.  (218)  into  Eq.  (205) 

W{x,  k,t)  —  —  [  S*(k  +  ±0,  0)  S(k  -  9/2, 0)  e-i9xeiHk+9/2)-u(k-e/2)]t  dQ  (231) 

2vr  J 


Starting  with 


S*  (, k  +  \e,  0)  S(k  -  .§0,  0)  e~idx  do 


W{x,k,  0)  ~  2k  !  ^  yv  '  2” ’ '1>  ^ v'1'  21 


we  invert  to  obtain, 


(232) 


S*(k  +  \e,  0)  S(k  -  \e,  0)  =  /  W(x,  k,  0)  ei9x  dx 


(233) 


and  inserting  into  Eq.  (231)  we  have 

W(x,  k,t)  =  ^  jj  W(x',  k,  0)  e-i°W-x)eimk+o/2)-w(k-e/2)]t  dQ  dx,  (234) 

This  is  exact.  It  is  convenient  to  write  this  in  terms  of  a  Green’s  function  for  the  Wigner 
distribution, 

i  r 

(235) 


where 


W(x,  k,  t)  —  —  I  W (. x ',  k,  0 )L(x’  —  x,  t )  dx' 
2vr  ./ 


L{x'  —  x,  k,  t)  =  /  e-i°W-x)eilw(k+e/2)-w(k-e/2)]td9 


(236) 


It  is  very  interesting  to  find  an  approximation  to  the  Wigner  distribution  with  the  follow¬ 
ing  aim.  The  evolution  may  be  simple  because  we  are  in  phase  space  and  it  is  generally  the 
case  that  evolution  equations  take  on  a  particularly  simple  form  in  comparison  to  evolution 
of  the  density  itself.  Expand  [W(k  +  9/2)  —  W(k  —  9/2)}  in  9 

pr/(2n+l)/i,\  f)2n+l  i 

W(k  +  e/2)-W(k-e/2)  =  YI  (2n+1)i  2„  ~v(k)t)  +  -vW(k)93  ...  (237) 


n=0 


where  if2n+1\k)  is  the  2 n  +  1  derivative  with  respect  to  k.  Keeping  only  the  first  term 


L(x'  -x,t)~  /  e-ie^-x)eiv{k)td  d9  =  5(x'  -x  +  v(k)t) 


(238) 


and  substituting  into  Eq.  (119)  we  have 


W(x,k,t)  ~  W(x  —  v(k)t,k,  0)  (239) 

This  is  a  remarkably  simple  result  and  gives  a  very  simple  way  of  propagating  the  Wigner 
distribution  for  pulses.  Note  that  no  calculations  are  required,  just  a  simple  substitution. 
That  is,  no  equations  have  to  be  solved  or  integrated.  This  method  may  be  a  good  practical 
method  for  evolving  pulses  in  time.  We  now  address  how  the  method  can  be  improved. 
What  is  needed  is  further  approximation  by  keeping  more  and  more  successive  terms  in  the 
expansion  given  by  Eq.  (123).  Going  to  the  next  term  we  have 


L(x'  -x,t)~  J  e-ie{x'-x)eiv(k)te+iv"{k)te^/2A  dQ 

(240) 

—  j  e-i0[(x'-x-v(k)t]eiv"(k)t63/24  dQ 

(241) 

This  integral  can  be  expressed  in  terms  of  Airy  functions  but  we  do  not  do  so  here.  This 
may  offer  a  considerable  better  approximation  than  Eq.  (239) 

We  point  out  that  the  recovery  of  the  pulse  from  the  Wigner  distribution  is  an  important 
problem.  The  signal  can  be  recovered  from  the  Wigner  distribution  up  to  a  constant  phase 
factor.  However  in  this  case  the  constant  phase  factor  may  be  a  function  of  time!  This  issue 
has  been  considered  by  Leavens  and  Mayato. 

Example.  Consider  the  initial  pulse  given  by 

u{x,  0)  =  (q;/7t)  1/4  e~ax2 /‘Z+iPx2 /2+ikox 

(242) 

The  Wigner  distribution  is 

W(x,k,0)  =  -  e-a*2-(fc-/J*-fco)7a 

7 r 

Therefore  we  have 

(243) 

W (x,  k,t)  ~  W(x  —  v(k)t,  k,  0) 

1  ^--Oi(x—v(k)t)2  —  (k—kQ—ft(x—v(k)t))2/a 

7T 

(244) 

(245) 

This  is  a  general  result  for  the  initial  pulse  where  the  dispersion  relation  is 

now  chose 

arbitrary.  If  we 

W(k)  =  ck  +  y/c2/2 
v(k)  =  c  +  y/c 

(246) 

(247) 
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we  have 


W(x,k,t)  ~  W(x  —  v(k)t,k,  0)  (248) 

=  W(x  —  (c  +  k'yt),  k,0)  (249) 

=  —  exp  I"— a((x  —  c  —  kjt)2  —  (k  —  j3{x  —  c  —  k-yt  —  k0)2/a]  (250) 

7 T 

Explicitly 

Wa(x,k,t)~  —  exp  [— a((x  —  c  —  k'yt)2  —  (k  —  /3(x  —  c  —  k'yt  —  k0)2/a]  (251) 

7 T 

It  is  remarkable  that  for  this  example  the  answer  is  exact! 

A  further  example  is 

W(k)  =  ck  +  7/c3/3  ;  v{k)=cJr^k2  (252) 


and  hence 


Wa(x,k,t )  =  W (x  —  v(k)t,  k,  0) 

=  VF(x  —  ct  —  7 k2t,  k,  0) 

^—a((x—ct—'yk2t)2—(k—ko—l3(x—ct—'yk2t)2/a 

7T 


(253) 

(254) 

(255) 


8  Case  B 


We  now  address  the  case  where  we  are  at  a  fixed  position  and  generate  a  pulse,  which  is  the 
situation  appropriate  in  radar,  sonar  and  fiber  optics,  etc.  We  generate  w(0,  t )  where  we  have 
taken  x  =  0  where  the  pulse  is  being  generated.  Many  new  concepts  are  now  introduced  for 
this  case.  An  important  issue  is  whether  we  stick  with  the  definitions  of  the  relation  of  pulse 
and  spectrum  as  before  or  we  change  to  make  it  easier  and  more  conventional  to  standard 
signal  analysis.  We  have  decided  to  take 


(256) 


for  the  fundamental  solution.  Substituting  into  the  wave  equation 

^  dnu  ;  dnu 

2.a^  =  2>^ 

n= 0  n= 0 

one  obtains 


N  M 

\n 


£' 


dnylM 
n= 0  n= 0 


-  J2h(-ik)n  =  0 


(257) 


(258) 
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Solve  for  k  in  terms  of  c o 


k  =  K{u) 

(259) 

The  general  solution 

is  then  given  by 

u(x,t )  =  _  f  F(u,x)  eluJt  du 

v27t  J 

(260) 

with 

F(uj,  x)  =  F{uj,  0)  e~iK{u})x 

(261) 

where  F(u>,  0)  is  the  spectrum  evaluated  at  x  —  0, 


F(u,  0) 


J  u(0,  t )  e  luJt  dtt 


(262) 


Also, 

F{ u,x)  =  _  f  u(x,t)  e~luJt  dt  (263) 

v27r  J 

That  is  u(x,t )  and  F(u,x)  form  Fourier  transform  pairs  for  any  x.  In  analogy  with  group 
velocity,  K'{u),  will  be  of  importance  .  In  particular  we  define  the  frequency  transit  t(u) 
for  a  given  frequency  by 

t(u)  =  K\u)  (264) 


We  call  it  the  frequency  transit  time  because  as  we  will  see  it  will  turn  out  to  be  the  time  taken 
for  a  given  frequency  to  travel  a  unit  distance.  We  now  discuss  some  important  concepts  for 
this  case  and  ones  thinking  must  change  from  the  previous  case.  One  must  visualize  that 
one  is  staying  put  at  a  particular  place  and  is  measuring  a  number  of  important  quantities. 
Suppose  we  have  a  signal  given  by  fit),  the  mean  time  is  traditionally  defined  by 


(t)  =  Jt  \f(t)\2dt  (265) 

The  only  difference  between  this  and  the  usual  definition  is  that  now  will  be  a  function  of 
position.  Similarly  the  duration  of  the  signal  is  given  by 

<?2  =  J (t  ~  (t))2\f(t)\2dt  (266) 


In  the  case  of  pulse  propagation  we  use  the  same  definitions  but  this  signal  will  be  u(x,t ) 
and  hence  is  a  function  of  position.  We  now  develop  the  analogous  ideas  and  formulas  for 
this  case  as  we  did  before  and  will  introduce  a  number  of  new  and  interesting  viewpoints. 
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We  first  developed  the  general  approach  in  terms  of  calculating  moments,  that  is  moments 
at  a  particular  point  in  space, 


(tn)x—  I  tn  \u(x,  t)\2  dt 

=  I  F*(u,x)TnF(u,x)du 

where  T  is  the  time  operator  in  the  frequency  domain 

d 


T  =  % 


du 


(267) 

(268) 


(269) 


We  now  discuss  each  moment  and  its  physical  interpretation.  Consider  the  first  moment.  It 
is  given  by 

( t  )x  =  ft  |  u(x,  t)  | 2  dt 


and  this  evaluates  to 


=  /  F* (u ,  x)T F (u ,  x)  du 


( t)x  =  {t)0  +  Tx 


(270) 

(271) 

(272) 


where 


T  —  I  r(u)\F(u,0)\2du  (273) 

These  equations  lead  to  the  following  interpretation.  r( u)  is  the  time  that  it  takes  a  sine 
wave  generated  with  frequency  u  to  travel  a  unit  distance  and  T  is  the  average  time  for  all 
frequencies  where  the  spectrum  is  used  as  the  weighting  function.  Therefore  the  mean  time 
at  position  x,  ( t)Xl  is  the  mean  time  at  position  (t) o  plus  the  time  it  took  to  travel  from 
the  initial  position  to  x.  Now  that  quantity  is  Tx  because  as  just  discussed  T  is  the  time 
per  unit  distance  and  x  is  the  distance  traveled  since  the  particle  started  at  x  =  0  initially. 
For  the  second  moment 


where 

Where 


(t  )x  =  (t/o  +  x(tT  +  Tt))0  +  x2(t2)o 


(tT  +  Tt))  o 


(rT  +  Tt))x  =  J  F*(u>,x)  [t(u)T  +  Tt(u)\  F(u,x)du 

f  d  d 

=  /  F*( u,x)  r(u>)v(k)i— — \-i—r(u>)v{k)  F(u,x)du 

—  if  F*(uj,x)  1 2  r(u) — ^  7  '  +  r'(u)F(k,  t)  \  du 


du 


(274) 

(275) 

(276) 

(277) 

(278) 
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A  standard  notation  is 


(279) 


(280) 

(281) 

(282) 

(283) 

(284) 

(285) 


then  in  general 


( Tt)x  =  \(tT  +  Tt)x 

(286) 

=  -J  r(u)^-\F(u,x)\2 du 

(287) 

Now  let  us  consider  how  duration  changes  with  position.  As  before  we  note  the  dominant 
term  is  x2a 2  (for  large  x )  and  its  coefficient  is  manifestly  positive.  Therefore  for  large 
distances  the  duration  must  go  to  infinity  no  matter  what  the  duration  is  at  the  point  where 
it  is  generated.  That  is  at  large  distances,  and  in  particular  at  infinity  the  duration  is  infinite. 

We  now  calculate  at  what  spatial  points  xc  the  duration  is  smaller  than  the  duration  at 
the  point  of  generation,  x  —  0.  The  duration  will  be  shorter  if  a^x  <  cr|0,  which  gives 

CovL0 

0  <  xc  <  -2 - 

°T 

(288) 

It  will  be  a  minimum  at  the  following  point 

Cov^io  1 

xm  —  9  —  n  XC 

oy  2 

(289) 

and  the  duration  at  %  is 

f)  o  COV+^-It*  r\  / 

(Jt\XM=°t  |0  a2  =ffi1o(1  PtT  1 

x)  (290) 
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where 


(291) 


2 

Ptr  |0 


CovtT|0 


|0 


The  Covariance  Between  Time  and  Transit  Time.  As  before  the  covariance  is  the  crucial 
quantity,  but  now  it  is  the  covariance  between  time  and  frequency  transit  time 


CovtT|x  =  {tr)x-  (r)x(t  )x 

If  we  express  the  spectrum  in  terms  of  amplitude  and  phase  then 

F(w,x)  =  |  F(u},x)\ei'Ku’x) 

=  F(u>,  0)  e-iK(uj)x 
=  |  F(w,0)|e^(w’o)e-<K(w)* 

and  therefore  the  phase  as  a  function  of  position  is 

u,  x )  =  0)  —  K{oj)x 

and  hence 

dip{u,x)  dip(u,  0)  , 

= —&r  ~  K  ^ 

or 

dtl>(u,x)  dip(u,0) 

—faj  =  —faj  TX 
Substituting  this  into  Eq.  (287)  we  have 

(t  t)x  =  (t  r)0  +  (  t2  )x 

Therefore 

Cov^i*  =  ( tr)x  -  {r)x(t)x  =  (tr)0  +  (t2)x  -  (r)0((t)0  +  {r)0x) 


(292) 

(293) 

(294) 

(295) 

(296) 

(297) 

(298) 

(299) 

(300) 


giving 


The  correlation  coefficient  is 


C'OVtrla;  CoV^t|q  T  X  CT 


Ptr\x 


CoVtr|a 


Covtr|o  +XCT 2 


and  we  have  that 


at\xaT\x  crr|0 yj cr^Q  +  2x  CovtT|0  +  x2al 

Ptr\x  *  1 


x  — >  oo 


(301) 


(302) 


(303) 
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Covariance  Between  Time  and  Frequency.  It  is  interesting  to  calculate  the  covariance  be¬ 
tween  time  and  frequency  at  a  given  position, 

Covto| x  =  (tuj)x-  {t)x{u)x  (304) 

Now 


(tu;)x  =  \{ uT  +  Tu)x 

(305) 

f  dif{u,x)  2 

=  —  /  u - - - \F(u),x)\  duo 

J  OU 

(306) 

—  —  ou[  — - - tx]\F(u,  x) \  du 

(307) 

f  jdip(u,  0)  2 

=  —  u  \  — - - tx\\T  (u,  0)  dco 

J  OOJ 

(308) 

Hence 

(tuj)x  =  (tu)o  +  (u  T  )0X 

(309) 

and  therefore 

Covtu\x  =  {ut  )x  -  (u)x{t)x 

(310) 

=  (tu) 0  +  (cvt)0x  -  (i 'jj)o((t)o  +  Tx ) 

(311) 

—  Covta,|0  +  x  CovWT|0 

(312) 

where  as  usual 

CoVWT|o  =  {u  T  )o  -  {u  )0(  T  )0 

(313) 

Relation  between  frequency  transit  time  and  group  velocity.  If  one  interprets  the  group  ve¬ 
locity  as  the  velocity  that  a  certain  frequency  has  and  if  one  interprets  the  frequency  transit 
time  as  the  time  it  takes  a  frequency  to  travel  a  unit  distance,  then  clearly  there  should  be 
a  relation  between  the  two.  Let  us  think  in  terms  of  the  analogy  with  runners.  Suppose  the 
velocity  of  a  runner  is  v  and  hence  in  a  time  T  he  travels 

x  =  vT  (314) 

Now  his  frequency  transit  time  r  is  the  time  it  takes  to  travel  a  unit  distance,  that  is  given 
by 

r  =  T/x  =  1/v  (315) 

and  hence  we  expect  the  relation  between  frequency  transit  time  and  group  velocity  to  be 

r(  u>)  =  l/v(K(u)).  (316) 

That  is  indeed  the  case. 
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8.1  Instantaneous  Frequency  and  Local  Spatial  Frequency 

As  before  we  define  amplitude  and  phase  by, 

u(x,t)  =  A(x,t)eitp^  (317) 

but  because  of  the  change  in  convention  for  this  part  we  define  instantaneous  frequency  as 

Ui(x,t)  =  (318) 

and  spatial  local  frequency  by 

ki(x,t)  =  ~^tp{x,t)  (319) 

Asymptotic  solution.  We  now  give  the  asymptotic  solution  for  case  B.  One  sets 

K'(us)  =  t/x  (320) 

and  hence 

ua(x,  t )  ~  F(u,  (321) 

and  the  amplitude  and  phase  are  given  by 


\ua(x,t)\  =  \F(»*0)\J  xK]{u)  (322) 

<pa(x,  t)  =  0)  —  K{uj)x  +  cot  —  7rsgn  K" /4  (323) 

Instantaneous  Frequency.  To  obtain  the  instantaneous  frequency  we  have 


d<pa(x,t) 

Wi(x,t)  = 

(324) 

duo  \dip  dK(u>)  1 

-  dt[du,  X  du  +tJ  +  " 

(325) 

which  gives 

.  .  I" dw  dip~\ 

""(x't)=  [aFU  j+" 

(326) 

and  since 

(327) 

giving 

du  1  K'(k) 

dt  ~  xK"(u)  tK"(k ) 

(328) 
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Hence, 


Ui(x,t)  = 


dip 


xK"{u )  du 
K'(k )  dip 


+  c 0 

+  u 


tK"(k ) 

Spatial  local  frequency.  The  spatial  instantaneous  frequency  is 


d 

ki(x,t )  =  -  —  (pa{X,t) 
dx 


dip  dus 


dus  dx 
dip 


+  t 


du. 


dus 

dip 


+ 1  —  x 


+  t  —  x- 
dus  x 


dx 
dK{us ] 
dus 
dujg 
dx 


K{us)  -  x 


dK{us )  dus 
du ,  <9a; 


<9a; 


+  K  (us) 


t 


+  A  (^s) 


or 

&iOM) 

But  from  Eq.  (320)  we  now  have  that 


dip  dus 
dus  dx 


+  K  ( us ) 


giving 


K"(us ) 


(9a: 


—t/x2 


dus  t 

dx  x2K"(us ) 


(329) 

(330) 

(331) 

(332) 

(333) 

(334) 

(335) 

(336) 

(337) 

(338) 


and  hence 


ki(x,t ) 


t  dip 
x2K"(us)  dus 


+  A  (us) 


(339) 


8.2  Wigner  Distribution  (Time-Frequency) 

The  Wigner  distribution  now  is  the  conventional  time-frequency  one.  Again  we  will  use  the 
standard  notation  for  the  signal  but  now  for  the  spectrum  we  have 


u(x,t)  =  A(x,  t)eiip(x,t^ 

(340) 

F(u,x )  =  _  f  u(x,t)  e~luJt  dt 

v27t  ,/ 

(341) 

=  B{  u,x)ei^u}’x) 

(342) 
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The  time-frequency  Wigner  distribution  is  then 


W(t,  u,x)  =  —  I  u*  (x,t  —  It)  u(x,t  +  \t)  e  lTUJ  dr 

Z7T 


—  — —  /  F*{uj  +  \6,x)F{u-\0,x)e-mde 


2vr 


The  instantaneous  frequency  is  the  first  conditional  moment  and  is  given  by 


(343) 

(344) 


Mt,*  =  ujW (t,  id,  x)du 

d  .  . 

=  at'p{f’x) 

and  the  group  delay  is 

(*>u,,®  =  B(u,x)  J  tW(t,u,x)dt 

Also  the  covariance  between  time  and  frequency  is 

(cut)x=  j  tcuW(t,u>,x)dtdcu 
d(p(t,  x 


=  t- 


=  -  UJ 


d t 
dip(u),  x 


du 


u(t,  x)\2dt 

WM\2 


Now  for  pulse  propagation  we  have 

F(cj,x)  =  F(cj,  0)  e~iK^)x 

=  B(u,  0)e*^(aJ,0)_iA  ^)3: 

and 


(345) 

(346) 

(347) 

(348) 

(349) 

(350) 

(351) 

(352) 

(353) 


which  gives 


B(u,  x)  =  B(u,  0) 

(354) 

■e- 

e 

■e- 

jT 

o 

1 

IT 

(355) 

d 

(356) 

=  0)  +  K'(u)x 

(357) 

d 

=  —Q^ij>(u},0)  +  t(uj)x 

(358) 
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At  x  =  0 

O 

1 

■e- 

1T 

O 

(359) 

and  therefore 

<*>o,,*  =  (*L,0  +  T(u)X 

(360) 

and  also 

(ut)x  =  J u  1  F{u,x)  |2 

(361) 

=  -/  u(-^ip(u,0)  -t(u)x)\F(u,0)  I2 

(362) 

=  {ut)o  +  (ut(u))0x 

(363) 

which  agree  with  the  results  of  Part  B. 

Now  consider  the  relation  of  the  Wigner  distribution  at  position  x  and  how  it  is  related 
at  position  zero.  We  write 


w(t,  lu,x)  =  —  /  f*(lu  +  \b,  o)  f(uj  -  e/2,  o)  de  (364) 

2tt  I 


and  also 


which  gives, 


W(t,  u,0)  =  —  /  F*(u  +  \B,  0)  F{ u  -  \9,  0)  e~m  d9 


F*  (u  +  \9,  0 )F{(j-  \9,  0)  =  /  W(t,  id,  0)  em  dt 


Inserting  this  into  Eq.  (364)  we  obtain 


which  as  before  we  write  as 


(365) 


(366) 


W{t,  lu,x)  =  —  W(t',  u,  0)  e-W-t)ei[k(u+9/2)-k(u,-e/2i }x  dB  dt,  (367) 


W(t,  u,  x)  =  —  W {tr,  lu,  0 )L(t'  —  t,  x)  dt ' 
2vr  ./ 


(368) 


Lit'  -t,Ld,x)  =  j  e-W-t)ei[K(,u,+0/2)-K(u,-0/2)]x  de 
Now  expand  [k(ou  +  9/ 2)  —  k( u  —  9/2)]  in  a  power  series  in  9  to  obtain 

00  jy(2n+l)  (.  .\  n2n+l  i 

K(u,  +  9/2)  -  K{*  -  9/2)  =  X  ~  T{*)e  +  -r<2V)93 

n= 0  '  ' 


(369) 


(370) 


where  K^2n+1^(ct/)  is  the  2 n  +  1  derivative  with  respect  to  to.  Keeping  only  the  first  term 


L{t'  -t,x)~  /  d9  =  5(t'  -  t  +  t(u)x) 


(371) 
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Substituting  this  into  Eq.  (368)  we  have 


W ( t ,  u,x)  ~  W{t-  t(lv)x,  u,  0) 

As  an  example  consider 

u(0,t)  =  (a/w)  V4e-ai2/2+m2/2+iuJ0t 
and  the  Wigner  distribution  is 

W(t,u,  0)  =  I  e-ati-(u-pt-u>0)i/* 

7 r 

and  therefore 

W (t,  U>,x)  ~  W(t  —  t(u>)x,  U),  0) 

^—a(t—T(uj)x)'2  —  (uJ—u>o—l3(t—T(u))x))2/a 
7 r 

If  we  further  assume  that 

K(u)  =  Duj  +  7o;2/2 
r(  uj)  =  D  +  70; 

We  have 


(372) 

(373) 

(374) 

(375) 

(376) 

(377) 

(378) 


W(t,u,x)  ~  W(t  —  t(u)x,  lu,  0)  (379) 

=  W(t  —  (D  +  urfx),  u,  0)  (380) 

=  —  exp  \—a(t  —  (D  +  uj^)x)2  —  (c 0  —  (3{t  —  D  —  u'yx  —  cao)2/al  (381) 

7T 

9  Space-Time  Signals  and  Distributions 

9.1  Time-Frequency/Spatial-Spatial  Frequency  Representations 

In  the  previous  sections  we  considered  a  joint  representation  for  time  and  frequency  for  a 
given  position.  We  also  considered  joint  space  and  spatial  frequency  distribution  for  a  given 
time.  However  as  of  now  there  is  no  way  to  handle  relations  that  involve  any  mixed  variables, 
for  example  position  and  frequency.  In  this  section  we  show  how  one  can  define  for  pulse 
propagation  a  four  dimensional  representation  involving  jointly  the  four  quantities,  time, 
frequency,  space,  and  spatial  frequency  and  for  which  the  previous  ones  are  special  cases. 
We  will  develop  the  material  in  this  section  in  the  following  way.  First  we  will  derive  the 
equations  keeping  everything  general  and  then  specialize  to  the  pulse  propagation  case. 
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9.2  General  results 

What  we  seek  is  [7] 

W (. x ,  k,  t,  u)  =  Joint  distribution  of  time,  frequency,  space  and  spatial-frequency  (382) 
We  use  the  four  dimensional  Wigner  distribution  defined  by 

W(x,k,t,u)=  jj  u*(x- \tx^~  ^T)u(x+^Tx,t  +  ^T)e~lTxk-1™  drdTx  (383) 

Now,  the  four  dimensional  Wigner  distribution  may  be  written  in  a  number  of  different 
ways  all  of  which  are  useful  depending  on  the  calculation  being  done.  First,  we  define  the 
two  dimensional  spectrum  of  the  signal,  G(k,  u) ,  and  also  list  our  previous  definitions 

u(x,t)  =  A^x^ty^  (384) 

S(k,  t )  =  — L  [  u(x,  t )  e~ikx  dx  =  B(k,  f)e^(M)  (385) 

v2vr  J 

F(u,  x)  =  — L  [  u(x,  t )  e~juJt  dt  =  B(u,  x)e^u’x)  (386) 

v27T  ,/ 

G(k,  y-lJJ  u{x,  t )  e~jut~jkx  dt  dx  =  L(k,  u)  (387) 

With  these  definitions  we  have 

W(x,k,t,(j)=  (J^j  II  G*{k  +  ^ex,u  +  ^d)G(k-^ex,u^^e)e-^x~^dffdex  (388) 
and 

W(x,k,t,u)=  (J^j  II  S*(k  +  ^9x,t-\r)S(k-^dx,t-\r)e-^x-^d9dT  (389) 
and  also 


W(x,k,t,u)=  jj  F*(u  +  ^9,  x  -  ^tx)F(cj  -  ^9,  x  +  ^rx)e  Jte  JTxk  d9  drx  (390) 

The  marginals  of  the  distribution  are  derived  and  given  in  the  Appendix. 

Moments.  There  are  of  course  numerous  moments  that  can  be  developed  however  we  will 
give  only  relevant  moments  and  give  them  in  the  order  which  gives  the  most  insight  into 
pulse  propagation.  The  calculations  are  given  in  the  appendix.  Consider  first  (co)Xtk,t,  the 
instantaneous  frequency  for  a  given  spatial  point  and  spatial  frequency  at  a  given  time, 

(v)x>k,t=  [  cuW(x,k,t,v)duj  (391) 
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1  1 

2n  2 i 


s  (x 


2  Tx^) 


ds{x  +  \rx,t )  ds*(x-\rx,t) 


dt 


dt 


S(X+  ~Tx,t) 


e~jT*k  dr r 


(392) 

In  the  Appendix  we  give  different  expressions  for  this  and  other  quantities.  Now  consider  the 
instantaneous  frequency  irrespective  of  spatial  frequency.  It  is  obtained  by  further  integrating 
Eq.  (392)  with  respect  to  the  spatial  variable, 


(u)x,t=  /  toW(x,  k,t,cu)  dudk 


I  S)x,k,tdk 


(393) 

(394) 


and  this  evaluates  to 

(u)x,t  =  A2(x,  t )  J^(x,  t)  (395) 

which  is  the  standard  result  for  instantaneous  frequency.  Now  consider  (k)Xttt0J.  Because  of 
the  symmetry  of  our  definitions  one  can  readily  write  expressions  for  {k)Xju,t  and  (k)Xjt- 
Position.  The  average  position  is  similarly  given  by 


(  X  )  k,t,u) 


J  xC(x,  k,  t,  u)  dx 


1  1 

27t  2  i 


dS(k,  t  +  \t) 
dk 


(396) 


dS*(k ,  t  —  |r) 


9TJLZS(M+§t) 


e~jTUJ  dr 

(397) 


We  now  average  further  to  obtain  the  average  position  for  given  wave  number, 


{x)k,t  = 

J  x  C(x,  k,t,u)  dx  du 

(398) 

= 

B2(k,t )  J^(M) 

(399) 

We  now  further  average  to  obtain 

(x)t  = 

j"  B2(k,t )  dk 

(400) 

Covariance.  We  now  consider  the  covariance.  In  the  appendix  we  show 

d <p(x,  t )  d<fr(x,  t ) 


(ku)x,t  =  A  L(x,t)- 


dt 


dx 


(401) 


Calculation  tool  and  physical  interpretation.  One  can  think  of  -§i<p{x,  t )  and  -§^<p{x,  t )  as  the 
frequency  and  spatial  frequency  in  the  x,  t  representation.  That  is  we  associate 

d 


OJ 


dt 


v(x,t) 


(402) 
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,  d 

k  -»•  Q^p[x,t) 


(403) 


Using  this  idea,  answers  can  generally  be  written  down  immediately  without  recourse  to 
calculations.  For  example,  suppose  we  want  (k) x>t  ,  then  we  can  immediately  write 


d 

{k)x,t  =  A2(x,t)—(p(x,t) 

The  same  viewpoint  can  be  taken  in  the  spectral  domain.  In  the  spectral  domain 

d 

t  -*• 

Suppose,  that  we  want  to  calculate  ( xt)k,w ,  we  simply  write 

dip(k,u>)  &ip(k,u) 


(xt)kjU>  =  B-(k,u)- 


(404) 

(405) 

(406) 

(407) 


dk  du j 

Using  this  method  we  have  the  following  results.  They  can  be  checked  independently  by 
direct  calculation.  The  local  covariance  of  frequency  and  spatial  frequency  is 

'  d(p(x,  t)  dip(x,  t)  d<p(x,t )  d <p(x,t) 


Co vXit(ku)  =  A  (x,t) 


dt 


dx 


dt 


dx 


(408) 


and 


Co  v(ku)  =  A2(x,t ) 


d(p(x,t)  d(p(x,t)  d(p(x,t )  d(p(x,t) 


dt 


dx 


dt 


dx 


dx  dt 


(409) 


Also 


(xt)k,u  =  B2(k,u ) 


d'd(k,co)  (hp(k,uj) 


du 


dk 


(410) 


9.3  Application  to  Pulse  Propagation 

The  above  results  are  general.  We  now  specialize  to  the  pulse  propagation  case.  Pulse 
propagation  is  imposed  when  we  take 


S(k,  t )  =  S(k,  0)  e~iW{k)t 

=  B(k,0)eiqp{kfi)-iW{k)t 


(411) 

(412) 
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or 


B(k,  t)  =  B(k,  0)  (413) 

il>(k,t)  =  il>(k,0)-W(k)t  (414) 


Consider 


d 


(4.,  =  |B(M)|2  \~ip(k,0)  +  W'(k)t 


6(u  -  W(k)) 

=  |B(fc,o)|2|wM  +  iy'(fc)iU(^-iy(*;)) 


(415) 

(416) 


This  is  an  interesting  result  and  can  be  interpreted  as  follows.  The  average  of  the  pulse 
at  a  certain  frequency,  spatial  frequency,  and  time,  varies  linearly  with  time  and  with  a 
velocity  equal  to  v{k)\  but  the  only  frequencies  allowed  are  those  given  by  u  —  W(k),  that  is 
those  that  satisfy  the  dispersion  relation.  Of  course  this  is  intuitively  obvious  but  we  have 
developed  the  mathematics  to  describe  the  situation.  If  we  average  over  all  frequencies,  then 


(■ x)k,t=  {x)ktUuJduj  =  \B(k,0)\~  {x)k0  +  v(k)t 


(417) 


This  is  exactly  the  result  obtained  in  previously  using  the  two  dimensional  Wigner  distrib¬ 
ution.  Also,  further  integration  over  spatial  frequency  gives 


{x)t=  J(x)kttduj=  J  \B(k,  0)|2  |  (x)k0  +  v(k)t 

=  ( x)o  +  (v(k))t 

=  (x)o  +  Vt 


du 


(418) 

(419) 

(420) 


which  is  Eq.  (51).  But  here  we  have  obtained  it  in  a  very  direct  manner.  Now  consider 

{x)t,u; 


(Xhu,  =  I  (X)ktt„  dk  =  j  \B(k,0)\2  |(x)fc0  +v(k)t |  8(u>  —  W(k))dk 


To  simplify  this  we  assume  that  the  solution  to  the  equation  u>  —  W[k) 
Then 


(x)t,u  = 


MM) 

W'(k) 


\B(k,0)\\=g{lv) 


0  is  A; 


(421) 

K(u>). 

(422) 


and  this  can  be  similarly  interpreted. 
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10  Examples 

10.1  Spread,  contraction  time,  etc.  for:  u(x,  0)  =  («/ e~ax 


In  this  example  all  quantities  can  be  solved  for  exactly  and  hence  offers  a  good 
verifying  the  results  we  have  obtained.  For  the  dispersion  relation  we  take 

W  ( k )  =  ck  +  7/c2/ 2 

and  the  initial  pulse  is  taken  to  be, 

u(x,  0)  =  (c0r )1/4e-a3:2/2+^2/2+^ 


We  define 


rj  —  a  —  i/3 


so  that 

u(x,0)  =  (a/n)1/4  e-vx2/2++iko3 
and  also,  for  convenience  we  define 


a 


a  = 


P'  = 


:«2+/52) 

p 


(a2  +  /32) 


The  initial  spectrum  is  calculated  as  follows, 

1 


S(k,  0)  = 


( a/n )174 

Vv 

(a/ 7r)1//4 


u(x,  0)  e  *fca;  dx 

(k  -  k0)2 


exp 


exp 


2r) 


a(k  —  k0)2  ,P(k  —  koy 

X  - 


2 (a2  +  (32)  2(a2  +  P2)  \ 


\Ja  —  ip 

_  (Q-/tQ  a' (k— fco)2/2— i/3'(k—  fco)2/2 

a/q;  —  Jp 

and  the  time  dependent  spectrum  is  therefore 
S(k,t)  =  S(k,  0)e~iW(k)t 


a 


M: 1/4 


Vv 

(a /it)1/4 
\Jol  —  ip 


exp 


exp 


i(ck  +  y/c2/2)f 


(k  -  k0 )2 
2  // 


2(a2  +  /52)  2(a2  +  /32) 


2 /2+if3x2 /2+ikoX 

case  for 

(423) 

(424) 

(425) 

(426) 

(427) 

(428) 

(429) 

(430) 

(431) 

(432) 

(433) 

(434) 

(435) 
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At  the  initial  time,  t  —  0,  the  mean  and  standard  deviations  of  x  and  k  are 


x)o  =  0 

(436) 

k)o  =  ko 

(437) 

2  1 

a*'°  "  2a 

(438) 

2  a2  +  /32 

afc'°  "  2a 

(439) 

Also  the  covariance  is 


Covxv  =  —  (440) 

The  moments  are  calculated  very  simply  from  the  time  dependent  spectrum.  The  average 
group  velocity  and  averaged  square  are 


( v )  =  c  4-  jk0 

( v  >  =  7  — ^ - (c  +  7&o) 

and  therefore  the  standard  deviation  of  group  velocity  is 

*5  =  (v2g)  ~  (v)2 

2  a2  +  /32 

=  7  - 

1  2a 


Hence 


{x)t  =  (c  +  7&0)  t 

{ x2  )t  =  [(1  +  P  7 1)2  +  7 2«2  ^2]  +  (c  +  7 M  V 

and  the  spread  is 


^2|f  =  ^2|o  [1  +  2/3- 7 1  +  72  (a2  +  /32)  f 2] 

2  ,  9,7/5  2  2«2+/^2 

=  <7s|o  +  2*— +  *  7 


2a 


2a 


1 


=  — [l  +  2/37t  +  72(a2  +  /32)t2j 
=  ^  [(1  +  /?7^)2  +72«2^2] 


(441) 

(442) 

(443) 

(444) 

(445) 

(446) 

(447) 

(448) 

(449) 

(450) 
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It  is  also  interesting  to  write  ax\t  in  terms  of  ax |0 


ax\ t  —  ax\0 


-  ^-2 
~  ax\0 


1  +  2/3y t  +  72 


4cr*|o 


+  /32  t2 


1  +  2/37*  +  7253V 

7io 


(451) 

(452) 


The  reason  this  point  is  important  is  that  it  shows  that  the  spread  in  time  is  not  a  linear 
function  of  the  original  spread.  So  for  example  if  we  make  the  initial  spread  very  small  the 
spread  in  time  will  relatively  increase.  This  may  be  counterintuitive  but  it  is  correct  and  is 
a  reflection  of  the  uncertainty  principle. 


Contraction  and  Spread.  The  pulse  will  contract  for  times 

2/3  ,  (3 


tc  < 


(453) 


7(a2  +  IP)  7(ct2  +  P 2) 

This  can  only  happen  if  either  7  or  /3  are  less  than  zero,  but  not  both.  The  minimum  width 
is 


2  2 

®x\tM  Tr|0 


a 


a2  a2  +  /32<T*10 


(454) 


Note  that  /3  determines  whether  the  pulse  will  contract  and  it  has  to  be  negative  for  con¬ 
traction.  Also  note  that  the  covariance  is  negative  for  negative  /3. 

One  can  also  obtain  the  value  of  / 3  that  will  maximize  the  time  of  contraction.  That  will 
be  the  case  when  a  =  (3 


1  1  _  1 

to  < - a  tM  - 

7/3  27/3  2 

and 

2  _  \  2 

ax\tM  2C’x\t° 

Thus,  the  pulse  cannot  get  any  narrower  than  half  of  its  original  width. 


(455) 

(456) 


10.1.1  Exact  Solution 


Our  approach  has  allowed  us  to  calculate  the  above  quantities  exactly  and  the  calculations 
have  all  been  done  simply  and  using  only  the  initial  pulse  and  spectrum.  That  is  the 
advantage  of  our  method.  It  is  now  interesting  to  verify  these  with  the  exact  solution  for 
the  pulse.  To  obtain  the  exact  solution  we  calculate, 


u(x, t ) 


wJs(k’t)e" 

1  (o/tt)1/*  f  r  [k  —  k0)2 
Vv  i6XPi  27? 


i{ck  +  ^k2 /2)t  +  ikx 


(457) 

(458) 
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This  integral  can  be  done  but  the  answer  turns  out  very  complicated.  We  have  found  that 
by  making  a  simple  transformation  first  one  gets  a  more  tractable  answer.  In  particular 
letting  k  — >  k  +  k0  we  have 


.  ,  1 
U[x,t)=  - 

V  27T 

(a/7r)1//4 

Vv  J 

[  exp 

'  £2 

—  - - i[c(/c  +  fc0)  +  7  (A;  +  k0)2/2}t  +  i{k  +  fc0)x  (459) 

_ 

1 

" 

(a/ 7r)1//4 

^  J 

[  exp 

— k2{- - h  iyt/2)  +  ik(x  —  ct  —  yk0t)  —  i(ck0  +  yk02/2)t  + 

2't/ 

(460) 

1 

~  v7^ 

(a/ 7r)1//4 

Vv  J 

[  exp 

■  /j.2 

- (1  +  irjyt)  +  */c(x  —  ct  —  7/7/)  +  i/co(^  —  ct )  —  iyk02t/2 

2rj 

(461) 

1 

~  v7^ 

(a/ 7r)1//4 

Vv  J 

[  exp 

'  ^2 

- (1  +  irjyt )  +  */c(x  —  ct  —  yk0t)  +  ik0(x  —  ct)  —  iyk02t/ 2 

2  7j 

(462) 


Carrying  out  the  integration  one  obtains 

rj  (x  —  ct  —  k0yt)2 


,  .  (a/n)1^ 

u(x,  t)  =  ,  exp 


a/1  +  iyrjt 


1  +  iyrjt 


+  ik0(x  —  ct)  —  iyklt/2 


Using  the  fact  that 


1  +  i'yrjt  =  1  +  7  f3t  +  iyat 


we  have  that 


(463) 

(464) 


1  +  iyrjt\2  =  (1  +  /3yt)2  +  7 2a2 1 2 

=  2a<rJt 


(465) 

(466) 


Therefore 


u(x,  t) 


(a/7 r)1/4 
a/1  +  i'yrjt 


exp 


Also, 


7]  (x  —  ct  —  k0yt)2 

2  2^J 


(1  +  7  f3t  —  i'yat)  +  i/c0(x 


ct)  —  iyk\t/2 

(467) 


77(1  —  i'yrjt)  =  (a  —  i/3)  (1  +  7  f3t  —  i'yat) 
=  a  —  i/3  —  iyt(a2  +  /32) 


(468) 

(469) 


and  hence 


u(x,  t)  = 


(a/7 r)1/4 
a/1  +  i'yrjt 


exp 


(x  —  ct  —  k^t)' 

Aaal\t 


[a  —  i/3  —  iyt(a 2  +  /32)}  +  i/co(x  —  ct)  —  iyklt/2 

(470) 
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This  allows  one  to  express  the  answer  in  terms  of  phase  and  amplitude 


ue(x,t)  = 


exp 


(x  -  ( x)tf 
4(7 l\t 


exp 


(27r(7Jt)1/4 

Xx-(x)t)2{(5  +  {a2  +  (32)1t} 


Aaal\t 


where 


The  amplitude  and  phase  are 


\ue{x,t)\  = 


■  exp 


Pe  = 


(27rCT.r|J1/4 

(x-(x  )t)2  {P  +  (a2  +  /32b t} 


5  =  -  arctan 

2  1  +  Pyt 


(x  -  ( x)ty 


+  iko(x  —  ct )  —  ik o  yt)2  —  id 

ayt 


4cr2,. 

IC|t 


4a<T*|t 


+  fc0(x  —  ct)  —  ky  yt/2  —  h 


^  )t)2  +  +  ^  +  Hx  _  rf)  _  7</2  _  5 


4a;CTz|t 


Aa°l\t 


(x-  (x  )t)2  {P  +  ( a 2  +  /?2)7t} 

4afTx|t 


+  /c0(x  -  (x)t)  +  kl^t/2  -  5 


=  (x)tYjt^vl\t  +  k0(x-  (z}t)  +  fco7*/2  -  5 

We  also  note  that  the  above  can  be  simplified  in  a  different  way  if  one  uses 

x  —  ct  =  x  —  (x)t  +  kojt 

(x  -  ct)2  =  (x  -  ( x  )t)2  +  2{x  -  { x  )t)k0^t  +  klrft2 


(471) 

(472) 

(473) 

(474) 

(475) 

(476) 

(477) 


(478) 

(479) 


10.1.2  Asymptotic  Solution 

It  is  also  interesting  to  obtain  the  classical  asymptotic  approximation  and  to  compare  to  the 
exact  result.  We  have  that 


W'{k)  =c  +  ~fk 

(480) 

W"(k)  =  7 

(481) 

Setting 

W'{k)  —  c  +  'yk  —  x/t 

(482) 

we  have  that 

,  x  —  ct 

ks  = - 

7t 

(483) 
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Now, 


W  ( k )  =  ck  +  ,yk2 /2 


(484) 


=  —  (x-ct)(x  +  ct) 


and  also  now  which  gives  that 


ksx  —  W{ks)t  = 


{x  —  ctf 


(485) 


(486) 


where 


The  asymptotic  solution  is 


ks  -  k0  =  —{x  -ct-  k0jt) 
it 

=  —Xx  -  ( x)t ) 
it 


(x)t  =  (c  +  ko^t 


ua(x,t)  ~  S(ks,0)  J^~e  47r/4exp  i  ^  ^ 


(487) 

(488) 


(489) 


(490) 


on  (a/7r)1/4  r  (X  ~  ct  -  kolt)2' 

S(M) = ~nr exp  [ — — J 

(a/n)1^  \  a(x  —  ct  —  ko'yt)2  ,(3(x  —  ct  —  kojt)' 

=  CXP  [  2 (a2+/32)72t2  “  1  2(a2  +  (32)l2t2 


and  therefore 


(491) 


(492) 


(a/V)1/4  FT  ^in/4  \  a(x  -  ct  -  k0-yt)2  .P(x  -  ct  -  ko^ft)2  (x  -  ct)2' 

u‘{x’ () ~ V T  exp [“  2(a w  - 'Whv+'— J 

(493) 


KOM)|2  = 


1  (a/ii)lt2 

li  Ja 2  +  /?2 


qQe  -  (rr  )t)2 
(a2  +  f32)l2t2 


(494) 


o/(a-  (^)t)2  (x-ct)2 

^  ~  _/3  27V  +  “^T“ 

Ax  —  (x  )t)2  (x  —  lx  )t)2  ,  /  \  \  1 

=  272f2  +  2yt  +fco(^-(a;)t)  +  -A:07t 

=  (x)t)2  +A:0(x-(x)t)  +  ^7t 


(495) 

(496) 


(497) 
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where  we  have  used 


(498) 

(499) 


x  —  ct  =  x  —  (x)t  +  kojt 

(x  -  ct)2  =  (x  -  ( X  )t )2  +  2(x-{x  )t)korrt  +  kffit2 


10.1.3  Comparison  of  Exact  with  Asymptotic  solution 

We  now  compare  the  asymptotic  solution  with  the  exact.  In  comparing  the  magnitude  of 
the  asymptotic  solution  with  the  exact  we  see  that  they  are  the  same  if  we  approximate  the 
conditional  standard  deviation  for  large  times  by  keeping  only  the  quadratic  terms 

ax\t-*oo  =  ^  t(«2  +  P2h2t2]  (500) 

Whether  this  is  a  general  feature  of  the  asymptote  solution  needs  further  investigation.  That 
is,  can  the  asymptotic  solution  be  obtained  in  the  following  way:,  write  the  parameters  of 
the  initial  pulse  in  terms  of  the  initial  position  and  initial  spread  and  then  substitute  for 
them  the  exact  moments  we  have  shown  how  to  calculate.  Is  then  the  asymptotic  solution 
the  one  that  is  obtained  when  only  the  quadratic  term  in  the  spread  is  kept? 

10.2  Example:  u(x,  0)  =  5(x) 

Suppose  we  take  an  impulse  at  x  =  0 


u(x,  0)  =  S(x) 


(501) 


The  initial  spectrum  is  given  by 


S(k,  0)  —  — == 

v  Z71 

We  consider  the  case  where  the  dispersion  relation  is  given  by 


W  (. k )  =  ck  +  7/c2/ 2 


(502) 


(503) 


and  hence 


S(k,t) 


_  -i(cfc+7fc2/2)t 

\Z2n 


(504) 
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u(x,  t )  = 


S(k,  t)eikxdk 


V2n 


—  _  /  eik(x-ct)-i'yk2t/2  jj 

2vr  / 


2tt 


2n  V  i^{t 


exp 


2m^t 


exp 


.  (x  —  ctf 
27 1 

(x  —  ct)2 
1  2 ~t 


The  exact  phase  is  therefore 


= 


(x  —  ay 
27 1 


and  the  exact  instantaneous  frequency  is 

.  .  d  .  . 


(x  —  ct)  (x  —  ct )“ 


7 1  27 f2 

(x  —  Xq  —  ct) 


27  f2 
(a;  —  ct)  (x  +  ct) 
2^2 


[2ct  +  (x  —  ct)] 


and  the  exact  spatial  instantaneous  frequency 


d 

ki{x,t)  =  —(p(x,t) 
(x  —  ct) 

7 1 


We  now  obtain  the  asymptotic  solution.  Using 

x  —  ct 

k*  = 


we  have 


7 1 

W (ks)  =  (x  -  ct)  (x  +  ct) 


ksx  —  W  ( ks)t  = 


(x  —  cty 
27 1 


Ua(x,t)  ~  S(ks,0)n 


tW"(ks 


aik3x—iW (ks)t-i'iTsgn  A'"/4 


V2tt\I  7 1 


^  ^ iks(x—ct)—'yk2t/2—i'Kry/& 


i2ii^t 


exp 


.  (x  —  cty 

1  2jt 


(505) 

(506) 

(507) 

(508) 

(509) 

(510) 

(511) 

(512) 

(513) 

(514) 

(515) 

(516) 

(517) 

(518) 

(519) 

(520) 

(521) 
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But 


<Pa(x,  t ) 


(x  —  ct)'2 
27 1 


(522) 


The  relation  between  exact  and  asymptotic  is 


<Pa(x,t)  =  c p(x,t ) 

(523) 

The  instantaneous  frequency  (asymptotic)  is  therefore 

(asymptotic)  =  u>i(x,t)  (exact) 

(524) 

Now  let  us  use  the  equation  from  the  text  to  derive  the  instantaneous  frequency, 

^i( 

,  x  chb  x 

x’t)  =  t*w"(k,)dk,  +  w{k‘) 

(525) 

We  note  that 

dtb 

dks  =  ° 

(526) 

and  therefore 

ca, 

i(x,t)  =  W(ks) 

(527) 

=  2^2^  ct)(x  +  ct) 

(528) 

Also 

,  .  ,  1  di(j  , 

ki{x’t)  =  tw«(k,)dk,  +  k- 

(529) 

=  ks 

(530) 

(x  —  ct) 
it 

(531) 

Which  agrees  with  the  above 

10.3  Example:  u(x,  0) 

II 

On 

i~r 

1 

b 

0 

The  reason  we  consider  this  example,  even  though  we  have  just  done  5(x)  will  become 

apparent 

u(x,  0)  —  S(x  —  Xq) 

(532) 

The  initial  spectrum  is  given  by 

S(k ,  0)  =  lP~ikX0 

(533) 

53 


and  hence 


gffc  t)  —  _ e-ikxo-i{ck+'yk2/2)t 

y/2l T 


u(x,  t )  = 


V2n 


2tt 


J  S(k,  t)elkxdk 

eik{x-xo-ct)-i'yk‘2t/2  ^ 


1  /2tT 

=  2tt\J  Tyt  6XP 


■  exp 


.  (x  —  x0  —  ct )s 
27 t 

.  (x  —  x0  —  ct)2 

1  2 ~t 


(534) 

(535) 

(536) 

(537) 

(538) 


Which  of  course  could  have  been  written  down  immediately  since  the  solution  is  translation 
invariant.  The  exact  phase  is  therefore 

(x  —  x0  —  ct)2 


<p(x,t)  = 


2yt 


and  we  can  write  down  the  answers  immediately  by  letting  x  — >  x  —  Xq 

.  .  d  .  . 

Vi[x,t)  =  —^<p[x,t) 

(x  —  x0  —  ct)(x  —  x0  +  ct) 


2  7t2 


and  the  exact  spatial  instantaneous  frequency  is 

d 

ki(x,t)  =  —<p(x,t) 

(x  —  x0  —  ct) 
7 1 

But  now  we  obtain  the  asymptotic  solution.  Using 

,  x  —  ct 
Ay  = 


7t 

W( Ay)  =  (x  -  ct)  (x  +  ct) 


ksx  —  W  ( Ay)t  = 


(a:  —  ct)2 
2yt 


«aOM)  ~  ^(fcs,  0)^ 


tW"(ks 


0iksx—iW (ks)t-iirsgn  K" /4 


^  —iksXQ  /  iks(x-ct)—^k‘it/2—m^/A 

V  7^ 

(x  —  ct)2  x  —  ct 

exp  1 — — - * - — x0 


t27T7t 


2yt 


7t 


(539) 

(540) 

(541) 

(542) 

(543) 

(544) 

(545) 

(546) 

(547) 

(548) 

(549) 
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But 


.  .  (x  —  ct )2  x  —  ct 

^  v  ’  ;  2  jt  7 t 

(x  —  ct  —  Xq)2  —  Xq2 
2 jt 


(550) 

(551) 


Thus  we  see  that  one  can  not  get  the  asymptotic  solution  form  the  previous  case  by  letting 
x  — >  x  —  Xq.  The  relation  between  exact  and  asymptotic  is 


<Pa(x,t)  =  C p(x,t )  - 


(552) 


The  instantaneous  frequency  (asymptotic)  is  therefore 


a Ji(x,t)  (asymptotic)  —  u>i(x,t)  (exact) 


(553) 


ki(x,t)  (asymptotic)  —  ki(x,t)  (exact) 


(554) 


Now  let  us  use 


to  derive  the  instantaneous  frequency.  We  note  that 


(555) 


=  - x0 


(556) 


and  therefore 


u>i(x,t)  =  — —  (— xq)  +  W(ks)  (asymptotic) 
t2  7 

1  .  w  ,  XXn 

= -  ct)(x  +  ct)  -  w 


(557) 

(558) 


and  therefore 


(x  —  ct)(x  +  ct)  =  (x  —  Xq  —  ct)(x  —  Xq  +  Ct)  —  Xq  +  2xXq 


(559) 


1  xx0 

Ui(X>  t)  =  2^2  -  Ct)(X  +  Ct)  -  ^2 

1  .  ..  .  x02 

=  2^^X~X°-  Ct^X  ~  Xo  +  Ct)  ~  2^ 


(560) 

(561) 


which  is  the  same  as  Eq.  (552).  Also, 
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ki(x,t) 


d'l/j 


tW"(k8 )  dks 


+  ks  — 


1  .  ,  x  —  ct 

=  -I  “®0  +  — J- 

71  7i 

x  —  x0  —  ct 
7 1 


(562) 

(563) 

(564) 


which  agrees  with  the  above. 
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1 


(565) 


10.4  Example:  u(x,  0)  =  elk°x 

Suppose  at  x  —  0  we  take 

u(x,  0)  =  — ^—elkox 
v  2tt 

The  initial  spectrum  is  given  by 

S(k,  0)  =  5(k  -  k0) 

We  consider  the  case  where  the  dispersion  relation  is  given  by 

W  ( k )  =  ck  +  7/c2/ 2 

The  time  dependent  spectrum  is 

S(k,  t)  =  S(k  -  k0)e~i{ck+^k2/2)t 


and  the  exact  answer  is 

u(x,t)  =  ~^=  f  S(k,  t)eikxdk 
y/2T[  J 

_  ^  —  i(cko+'ykn/2)t+ikox 

The  exact  phase  is  therefore 

p(x,  t )  =  ~(ck0  +  rykl/2)t  +  k0x 


and  the  exact  instantaneous  frequency  is 

.  .  d  .  . 

Vi(z,t)  =  -—<p(x,t) 

=  ck0  +  7&o/2 
=  W(k0)  (exact) 

and  the  exact  spatial  instantaneous  frequency  is 

ki(x,t)  —  k0  (exact) 


(566) 

(567) 

(568) 

(569) 

(570) 


(571) 


(572) 

(573) 

(574) 

(575) 


We  now  obtain  the  asymptotic  solution.  From  before  we  have  that 

x  —  ct 


kv  = 


ksx  —  W(ks)t  = 


7 1 

(x  —  ct)2 
27 1 


(576) 

(577) 
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and  therefore 


ua(x,t)  ~  S(ks,  0), 


tW"(ks 


=  -  ko)  J- \  exp 

Vi  V  7 1 

=  7lS(X~ir  “ 

=  -^=<5(x  —  ct  —  ytk0)\/yiexp 
V* 


0iksx—iW  (ks)t— 27rsgn  K"  j  4 


.  (x  —  ct)‘ 


2yt 
(x  —  ct)2 
1  2yt 
.  (x  —  ctf 
1  2yt 


and  therefore  the  asymptotic  phase  is 


<Pa{x,  t )  = 


(x  —  ct )“ 
2yt 


=  kljt/2 


(578) 

(579) 

(580) 

(581) 

(582) 

(583) 


The  instantaneous  frequency  (asymptotic)  is  (note  that  we  must  differentiate  Eq.  (582)  and 
not  Eq.  (583) 


Ui(x,t) 


( x  —  ct)  (x  —  ct)2 
yt  2  yt2 

(x  —  ct)  (x  —  ct)2 
yt  2  yt2 

2(x  —  ct)ct  +  (x  —  ct)2 


2  yt2 

(x  —  ct)  (x  +  ct) 
2  yt2 

=  W(ks) 

=  W(ko) 


(i asymptotic ) 


(584) 

(585) 

(586) 

(587) 

(588) 

(589) 


and  for  this  case  we  have  that  the  exact  equals  the  asymptotic.  Also, 

,  ,  .  (x  —  ct)  .  .  . 

kAx,t)  =  c -  ( asymptotic ) 

yt 


=  cks  ( asymptotic ) 

=  cks 


(590) 

(591) 

(592) 


Now  we  use 


and  using 


Vi(x,t) 


x  dip 
t2W"{ks)dks 


+  W(ks) 


dip 

dks 


(593) 
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(594) 


we  have 


u)i(x,t )  =  W(ks)  ( asymptotic ) 

=  W(k0) 


Which  is  the  same  as  Eq.  (574). 

Also 

h(x’t)  =  m*miir,+k'  {asymptoUc) 

=  ks 
=  k0 


(595) 

(596) 


(597) 

(598) 

(599) 
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10.5  Example:  u(x,  0)  =  elPx2/2+zkox 

At  t  =  0  we  take 

u(x,0)  =  eiP*2/2+ikox 

The  initial  spectrum  is  given  by 


/  * 

.{k  -  /c0)2 

/rxp 

— - — 

L  2/3  J 

We  consider  the  case  where  the  dispersion  relation  is  given  by 

W(k)  =  ck  +  ^k2/ 2 

and  hence 


S(k,t)  =  \l~jzx  P 


.  (k  —  ko )2  .  , 

-c  2/3° '  -  iW(k)t 


(600) 


(601) 


(602) 


(603) 


u(x,t)  =  -j=j  S(k,  t)eikx 


1  1 


^2n  V~ZP 
1  1 


V2n  y/—i/3 
1  1 

yf—W 


1  1 


1  1 


\^2rr 

1 


exp 

exp 

exp 

exp 

exp 


i — — — - i(ck  +  ^k2 /2)t  +  ikx 

^P 

k2 

-i— ~ - i[c(k  +  ko)  +  7  (k  +  k0)2/2]t  +  i(k  +  k0)x 

^P 


(604) 

(605) 

(606) 


-k2(~W  +  Ht/2)  +  ik(x  -  ct  -  7 k0t )  -  i(ck0  +  ^k2/2)t  +  ik0x 


(607) 


k 2 

—  (1  +  iPjt)  +  ik(x  —  ct  —  'ykot)  +  ik0(x  —  ct)  —  i^k02t/2 
Zfj 


(608) 


k2 

— i—  (1  +  Pyt)  +  ik(x  —  ct  —  7  k0t)  +  ik0(x  —  ct)  —  i^ko2t/2 

2p 


:  exp 


a/1  +  7 pt  [  2  1  +  7 pt 

The  exact  phase  is  therefore 


.gOr-rt-y)*  +  _  rf)  _  i<t/2 


(609) 

(610) 


P  (x  ct  ko'yt)2  2 

<p(a:,  t)  =  - - 1+  ^ - +  W  -  ct)  -  7V/2 


(611) 
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and  the  exact  instantaneous  frequency  is 


,  d 

Vi(x,t)  =  -—<p(x,t) 


P(x  -  ct  -  k0jt)  p2(x  -  ct-  ko^tf  2 

=  (c  +  fco7) - 7^—^ - +  7 - o/1  ,  ..^9 - +  k°c  +  7*o/2 


(1  +  7/ft) 

(3{x  —  ct  —  ko'-ft) 


1  +  7  /ft 


(c  +  fc07)  +  7 


2(l+7/ft)2 
/3(x  —  ct  —  k0jt) 


2(1  + 7/ft) 


+  /cqc  +  7/cq/2 


and  the  exact  spatial  instantaneous  frequency  is 


(3{x-ct-kolt) 

l[,)  (1  +  7 /ft)  +  0 

Thus  it  still  remains  a  chirp  but  with  a  different  chirp  rate. 
We  now  obtain  the  asymptotic  solution.  We  have  that 

x  —  ct 

kc  = 


7 1 


(612) 

(613) 

(614) 


(615) 


(616) 


and  therefore 


W(ks)  = 
ks-k0  = 
ksx  —  W  ( ks)t  = 


(x  —  ct)  (x  +  ct) 

2  yt2 

x  —  ct  —  ko'ft 


7 1 

(. x  —  ct)2 
2  jt 


ua(x,t)  ~  S(ks,0), 


tW"(ks 


0iksx—iW  (ks)t— 27rsgn  K"  j  4 


where  we  have  used  the  fact  that 


e-nr/4 


1 


(617) 

(618) 

(619) 

(620) 

(621) 

(622) 

(623) 
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and  therefore  the  asymptotic  phase  is 


/ 3yt{x  —  ct )2  —  (x  —  ct  —  yt/co)2 
2/3rj2t2 

(. x  —  ct)2([3^t  —  1)  +  2(x  —  ct)jtko  —  (itk0)2 
2/3ry2t2 

(x-ct-^tko)2  (x  -  ct  -  'ytko)2  7t/c02  ,  .  . 

= - wP - + - - +  -  +  (x-ct-1tk0)k0 

/  x  x  dib  TTT/,  . 

“i{x’t)  =  mv^k7)Ik,  +  w{k-) 

dv  (k  —  /.'()) 


We  note  that 


d/cs  f3 

x  —  ct  —  fc07t 
Pit 


and  therefore 


Ui(x,t)  = 


x  x  —  ct  —  k0"ft  (x  —  ct)  (x  +  ct) 


=  —  x- 


7t2  P^t  2yf2 

x  —  ct  —  fc07t  —  ct)  ( x  +  ct) 


/?7  2t3 


Also, 


&iOM)  = 


+ 


dip 


2  7t2 


(i asymptotic ) 


tW"(A;s)  dA;, 


+  ks  ( asymptotic ) 


x  —  ct  —  k0^t  x  —  ct 


/?72t2  '  7t 

x  —  ct  —  kojt  x  —  ct  —  kojt 


P'jH2 


+ 


7t 


+  A7 


(624) 

(625) 

(626) 

(627) 

(628) 

(629) 

(630) 

(631) 

(632) 

(633) 

(634) 

(635) 
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(636) 


10.6  Example:  u(0,t)  =  eluJot 

We  consider  the  case  where  the  dispersion  relation  is  given  by 

K(u)  =  7a;2 /2 


and  at  x  =  0 

u(0,t)  =eiuJot 

(637) 

Its  spectrum  is 

F( u>,  0)  =  V2ttS(ou  —  cn0) 

(638) 

and  hence, 

F{u,  x )  =  y/2 nd(cu  -  ou0)e~ilu}2x/2 

(639) 

and  therefore 

u(x ,  t)  —  1  F(u,x)elu}tdu 

v2vr  ,/ 

(640) 

—  e-i,y^ox/2+ju>ot 

(641) 

and  we  see 

that  the  phase  is  given  by 

<f(x,  t )  =  —^(jJqX/2  +  Uot 

(642) 

which  gives 

.  d  .  . 

Ui{X,t)  =  -r£P\?,t) 

(643) 

=  LV0 

(644) 

and  the  spatial  local  frequency  by 

d 

h(x,t )  =  ^—(p(x,t) 

(645) 

=  7^o /2 

(646) 

We  now  use  the  equations  derived  for  the  asymptotic  answer. 

We  have  to  first  obtain  u>s 

K'(us)  =  t/x  =  7 ujs 

(647) 

and  hence 

t 

UJS  =  — 

7 X 

(648) 

Further  for 

our  case  we  have  that 

^  =  0 

000 

(649) 
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Now 


and  also 


Vi(x,t) 


1  dip 

xK"(u> )  du  UJ' 
t 

'fX 


(650) 

(651) 


ki(x,t)  = 


and  we  see  the  answers  do  not  agree, 
the  asymptotic  signal  itself 


x2K"(ujs  )  dujs  ^  ^ 

K{us) 

7^/2 

7 1 

2x 

We  now  explore  why  this  should  be  case. 


(652) 

(653) 

(654) 

(655) 
Consider 


ua(x,t)  ~  F(uj,0)< 
Let  us  hrst  work  out  that 


xK"(u) 


Oiojt— iK(uj)x— 27rsgn  K"  /  4 1 


\(jJ=UJS 


cut  —  K(u)x  = 


2^x 


and  here  the  asymptotic  solution  is 


(656) 

(657) 


ua(x,t )  ~  F(u>,  0)< 


xK"(u) 


Diiot— iK(uj)x— 27rsgn  K"  /  4 1 


=  V2i r§( - Wo)\  — exp 

'yx  V  X1 

Hence  because  of  the  delta  function  we  now  obtain 

.  .  t 

Vi{X,  tj  UJ q 

'fX 


2^x 


and  also 


Hx,t)  =  ^ 

72cn0 


(658) 

(659) 

(660) 

(661) 

(662) 


which  still  does  not  agree  with  the  correct  answer.  However  if  we  take  the  asymptotic  signal 


and  rewrite  it  as 


u 


a 


(. x ,  t )  ~  V2nS( 


t  [T 

r  t2  l 

- Uo)\  — exp 

i - 

7 x  V  X1 

2^x 

t  [T 

'  t2  ' 

- Uo)\  — exp 

'yx  V  X1 

2^x 

(663) 

(664) 
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then  one  would  get  the  right  answer  if  one  differentiates  the  phase.  What  this  example  shows 
is  that  a  further  clarification  is  needed  when  one  deals  with  delta  functions.  In  particular 
for  this  case  it  is  not  clear  which  should  be  done  first,  the  delta  function  substitution  or  the 
differentiations. 
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10.7  Example:  u(0,t)  =S(t  —  to) 

Suppose  we  take  an  impulse  at  x  —  0 


The  initial  spectrum  is  given  by 


w(0,  t)  —  i 5(t  —  to) 


F(u,  0)  =  -=e~™° 
\J  In 


and  hence 


the  exact  answer  is 


F(u,x)  = 

V  2n 


1 


u(x,  t)  =  /_  J  F(lu,  x)ezutdu 

e-i'yu}2x/2-iu}(t0-t)  ^ 


\/2n 

—  f 

2n 


1  /  2tt 

=  7T\- — exp 
2n  V 


/ 

2'yx 


and  we  see  that  the  phase  is  given  by 


V{x,t)  = 


(t  -  t0)‘ 
2^x 


which  gives 


.  d 

Ui{x,t)  = 

t-tn 


'yx 


and  the  spatial  local  frequency  by 


d 

ki(x,t )  =  -—(p{x,t) 

(t  -  tp)2 
2'yx 2 


We  now  use  the  formulas  derived  in  the  text 


K'{us)  =  tlx  =  juis 


(665) 

(666) 

(667) 

(668) 

(669) 

(670) 

(671) 

(672) 

(673) 

(674) 

(675) 

(676) 
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and  hence 


uy  = 


")X 


KM  =  Ms/2  =  2^2 


Now 


and  hence 


di(j(tu,  0 

du 


—  ~tn 


.  .  1  dip 

Mdx’t)  =  7K^r)^;+h,‘ 

—  — ( — ^o)  H 

x'y  7 x 

= - -  ( asymptotic ) 

'fX 


which  is  the  same  as  the  exact. 


.  >  T~~  +  KM 


x2K"(ojs)  duis 


t 


X2r) 


(—to)  + 


i 2 


2'yx2 


t 2  —  tto 

2'fx2 


consider  the  asymptotic  signal  itself. 


ua(x,t)  ~  F(u,  0)„ 


xK"(u) 


Diujt— iK(uj)x— msgn  K"  /  4 1 


Let  us  first  work  out  that 


ut  —  K(u)x  = 


2^/x 


and  here  the  asymptotic  solution  is 


ua(x,t)  ~  F(u,  0)< 


xK"(u ) 


Diujt— iK(uj)x— msgn  K"  /  4 1 


|a;=a;s 


l  /~T~  t2  t 


(677) 

(678) 

(679) 

(680) 
(681) 
(682) 

(683) 

(684) 

(685) 

(686) 

(687) 

(688) 

(689) 
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10.8  Example:  u(0,t)  =  {a/ 7r)1/4e  at‘2/2+jUJot 


For  the  signal 

u(0,t)  =  (a/7r)1/4  e~Qt2/2+juJ0t 


the  initial  special  spectrum  is 


(a/7 r)1/4 

F(  u)  =  — —f= —  exp 


a 


and  the  exact  solution  is 


u(x,  t )  = 


a 


/t t)1/4 


V2a 

The  phase  and  amplitude  are 


|a  —  77a: 


exp 


(cn  -  cap)^ 
2a 


_^o  +  (^  +  7f) 


2a  4(^-  -  772:) 


,  .  —jxt2—uj0t/a2  +  uJo'yx/a2  1  27a; 

(p[x,t)  = - p— : - p - 1 —  arctan - 

4  [(2^)  +  72a:2]  2  a 


-7  a;a2f2  —  cnot  +  cno7a:  1  2^x 

- - - ^  9  „  „ - h  -  arctan - 

1  +  4a272a;2  2  a 


KM)  I  = 


a 


/7T)1/4  / 


V2 a  V«/4  +  72a;2 


exp 


1  /t2  —  4cuo7a;(t  —  ujq^x)  \ 

1  +  4a272a;2  / 


This  gives 


Ui(x,t)  = 


u0  +  2o?^xt 


1  +  4a272a;2 

This  is  a  chirp  even  though  a  pure  sine  wave  is  being  generated  at  x  —  0.  In  fact, 
uj  =  0  we  have  a  chirp. 


(690) 

(691) 

(692) 

(693) 

(694) 

(695) 

(696) 
even  for 


10.9  Example:  u(0,t)  =  e*^2/2+iw° l 

We  take 

u(0,t)  =  ei/3t2/2+iu,ot 

which  is  a  chirp.  We  have 


/  * 

.(w  -w0)2 

pew 

[  *  2/3  J 

and  hence, 


and  therefore 


F(w,x)  =  \  j-  exp 


.(u-UJof  .  2  ' 

-* - 2p - l^u  x /2 


u(x,  t )  = 


■ 


F(u;,:c)e^dw 

.(w  -  w0)2 
1  2/3 

V 


/  2tt/3  J 

f  exp 

I^J 

f  exp 

IWJ 

f  exp 

i^u2x/2  +  iut 


dw 


-  *7^  +  u0)2x/2  +  i(u  +  w0)t 


du 


.u1- 


exp 


1  +  P'yx 

and  we  see  that  the  phase  is  given  by 

Pit-  IUqx) 2 


-i— (1  +  P^/x)  +  iu(t  —  7  u0x)  —  i^u^x/2  +  iu0t 
Z/J 

'  ,p(t  —  ^U0x)2  .  2  ,  ' 

8  2(1  +  M  -  lw/2  +  I“°‘ 


du 


which  gives 


yfet)=~2(l+'^xi  ^W2  +  “' * 


9 


WiOM)  =  ^POM) 


t-^U0X 

=  '3TT7^  +  k'° 


and  the  spatial  local  frequency  by 


0 


fc*0M)  =  -^</?0M) 


P(t-~/U0x)  p(t-^U0xy  2 


7u;0  (1  +  /dyx)  +  (1  +  //bpr) 


7^5/2 


(697) 

(698) 

(699) 

(700) 

(701) 

(702) 

(703) 

(704) 

(705) 

(706) 

(707) 

(708) 

(709) 
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We  now  compare  to  the  asymptotic  solution.  We  have 


K'(u)  =  yu;  =  — 

X 

(710) 

and  that  gives 

t 

u  =  — 

JX 

(711) 

Therefore, 

F(u  =  — 

JX’ 

x  Cd/vr)174  r 

l=  VJ3  eXP[' 

( t  -  1Uqx)2~ 

7  2x2(3 

(712) 

and  hence 

.{t  -  'yu0x)2 
exp  i 

2y2arp 

t2 

+  i - i7rsgny/4 

z'yx 

(713) 

This  gives  an  instantaneous  frequency 


t 

UJi  =  — 


(714) 


10.10  Example:  u(0,t)  =  (a/7r)1^4e  atV2+*^i2/2+*woi 


This  example  is  mathematically  identical  to  the  example  previously  considered  and  hence 
we  do  not  give  the  details  but  just  the  results.  We  take 


K(u )  =  Du  +  7cu2/2 

(715) 

and  the  initial  pulse 

is  taken  to  be, 

u(0,t)  = 

(716) 

Defining 

rj  —  a  —  i/3 

(717) 

so  that 

u(0,t)  =  (a  /  7r) 1/4  e~vt2 /2++iu>ot 

(718) 

and  where  as  before 

,  a 

(a2  +  /32) 

(719) 

(3’  =  ^ 

(a2+/32) 

(720) 
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The  initial  spectrum  is 


and  further 


F(cu,0) 


(a /it)1/4 

Vv 

( a/n )1/4 
\Ja  —  if3 


■u(0, t ) 


exp 


exp 


e~iuJt  dt 

(u  -  cup)2' 
2r/ 

a(cu  —  cu0)2 
_  2(a2  +  /32) 


,/3(u;  -  cup)2 
?2(a2  +/?2) 


{ol/tt)1/4  e_Q'(a,_a)0)2/2-j/3'(w_a,0)2/2 
\/a  —  i(3 


F{u,x)  =  F(u,  0) 


(a/vr)1/4 

Vv 

(a/i t)1/4 
\Ja  —  i/3 


exp 


exp 


(cu  -  (Up)2  _ 
2  T) 

a(u)  —  cu0)2 
"  2(a2  +  (32) 


i(Du  +  7a;2 /2)x 


.P(u-u  o)2 

1  2  (a2  +  f32) 


i(Du  +  'yu2/2)x 


At  the  initial  position,  x  —  0,  the  mean  and  standard  deviations  of  t  and  lu  are 


(721) 

(722) 

(723) 

(724) 

(725) 

(726) 

(727) 


(t)  o  =  0 

(728) 

(cu  )o  =  cu0 

(729) 

2  1 

^  -  2^ 

2  a2  +  /32 

~  2a 

(730) 

(731) 

and  also 

Covfr  =  2i 

2a 

(732) 

The  average  transit  time  and  averaged  square  are 

(r)  =  D  +  7CU0 

(733) 

<t2}=72“22+/?2  (B  +  TU-o)2 

(734) 

and  therefore  the  standard  deviation  of  transit  time  is 

a2=(r2)-(r)2 

(735) 

2a2  +  (32 

1  2a 

(736) 
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Hence 


(t)x  =  {D  +  7CU0)  x 

( t2  )x  =  [l  +  2(3  7  x  +  72(a2  +  /32)  x2]  +  (D  +  7W0)2  x2 


(737) 

(738) 


and  the  spread  is 


Also, 


t2 

h|x 


=  °t\0  [l  +  2(3  y  x  +  y2(a2  +  (32)  x2] 

(739) 

2  ,  o  7/5  ,  2  2«2  +  /52 
=  at|o  +  2x2a+X7  2ffl 

(740) 

=  ^  [l  +  2(3  yx  +  72(ck2  +  (3 2)  x2] 

(741) 

=  7^  [(l  +  /37x)2  +  72a2a:2] 

(742) 

f2o  1  +  2Plx  +  72  f  4(t4  +  P2  j  2:2 

(743) 

f20  1  +  2(3yx  +  72  21  2:2 

(744) 

h|o 


Contraction  and  Spread.  The  pulse  will  contract  for  positions 

2(5  (5 


xc  < 


%  — 


(745) 


7<a2  +  ff2)  "  7<q2  +  ft2) 

This  can  only  happen  if  either  7  or  (3  are  less  than  zero,  but  not  both.  The  minimum  width 
is 


2  _  2  _COVMO 

at\xM  °*|0 


a 


a2  a2  +  (52<Tt^0 


(746) 


One  can  also  obtain  the  value  of  (5  that  will  maximize  the  time  of  contraction.  That  will 
be  the  case  when  a  =  (5 


xc  < - 7 

7  P 


Xm  — 


W  ”  2Ic 


and 


2  _  _  2 

®t\xM  7(^t\xo 


Exact  Solution.  To  obtain  the  exact  solution  we  calculate, 
1 


u(t,  x)  = 


F(u,  x)e 


iut 


1  (a/ 7r)1//4 

y/rj 


exp 


(a;  - 
2r] 


—  i(Du  +  yu2 /2)x  +  iut 


(747) 

(748) 

(749) 

(750) 
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and  hence 


u(t,  x ) 


1 

(cK  /  7T  ) 174 

v7^ 

1 

(oi  /  7T  ) 174 

y/2 7T 

1 

(cK  /  7T  ) 174 

1 

(a/ 7r)174 

v7^ 

V3 

exp 

’  ^2 

- - i[.D(c</  +  c^o)  +  c^o)2/2]x  +  +  c<Jo)£  (751) 

.  ^ 

exp 

x 

— cn2( - b  i'yx/2 )  +  iuit  —  -Dx  —  ^uj^x)  —  i(Duj0  +  7ca02/2)x  +  icu0t 

2  rj 

(752) 

exp 

"  O’2 

- (1  +  irj'yx)  +  7<x(t  —  -Dx  —  7ca0x)  +  7ca0(t  —  -Dx)  —  i^UQ2x/2 

2rj 

(753) 

exp 

"  O’2 

- (1  +  irj'yx)  +  iui(t  —  -Dx  —  7c a0x)  +  7cn0(t  —  -Dx)  —  i^UQ2x/2 

2  77 

(754) 


Carrying  out  the  integration  one  obtains 

77  (t  —  Dx  —  Uo'-yx)2 


(0/77) 174 

u{t,x)  =  ,  =  exp 

a/1  +  iyrjx 


Using 


we  have  that 


2  1  +  i'jrjx 

1  +  iyr/x  =  1  +  7  /fa  +  ijax 


+  ica0(t  —  Dx)  —  i^Uqx/2 


|1  +  iyr/x\2  =  (1  +  P'jx)2  +  rj2a2  x2 


=  2aa 


t\x 


(755) 


(756) 


(757) 

(758) 


Therefore 

(ck/tt)1/4 

u{t,x)  =  -7====  exp 


y/1  +  7770: 


77  (/  —  Da:  —  cao7x/ 
'2  2^1 


(1  +  7/fa  —  i'yax)  +  7<Xo(t  —  -Dx)  —  i^ulx/2 

(759 


Also, 


77(1  —  i'yqx)  =  (ck  —  i/3)  (1  +  7/fa  —  ijax) 
=  a  —  i/3  —  i^x(oi2  +  (32) 


(760) 

(761) 


and  hence 

.  .  (a/7r)1//4 

u{t,x)  =  .  exp 

Y  1  T  77772; 


(f  —  -Dx  —  ca07x)s 


[a  —  i/3  —  i'yx^a2  +  /3  )]  +  7cn0(t  —  -Dx)  —  i'jujQx/2 

(762) 
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and  also, 


ue(t,x)  = 


exp 


(* -(*)*)* 
4at\x 


exp 


(27r^|x)1/4 

.{t-  (t  )x)2  {P  +  (a2  +  P2)ix} 


4aat\x 


+  iu0(t  —  Dx )  —  icug  7^/2  —  i<5 


where 


The  amplitude  and  phase  are 


\ue{t,x)  |  = 


exp 


4(T|x 


<Pe  = 


(27rCTt|x)1/4 

(*-  (Ox)2{/5  +  (a2  +  /?2)7^} 


1  a^x 

o  —  —  arctan - - — 

2  1  +  p^x 


C t-(t)xY 


4a(7t\x 


+  cao(t  —  Dx)  —  eng  7^/2  —  <5 


+  («  -  (JM  (f  +  /32)7^  +  Mt  _Dx)  _  ^  ix/2  _  g 


4a(Tt\x 


4a(Jt\x 


(t-  {t  )x )2  {P  +  ( a 2  +  /32)7o;} 


+  ca0(t  -  (t)x)  +  <^07^/2  -  5 


=  -  ( t  >.r)2^  In  cr^  +  u;0(t  -(*}*)  +  ^o7^/2  -  <5 


(763) 

(764) 

(765) 

(766) 

(767) 

(768) 

(769) 


We  also  note  that  the  above  can  be  simplified  in  different  ways  if  one  uses 


t  —  Dx  =  t  —  (t)x  +  ca07a; 

(t  -  Dx)2  =  (t  -  ( t )  x)2  +  2  (t-  {t  )x)u0'yx  +  uffix2 

Asymptotic  Solution.  We  have  that 

W'(  u)  =  D  +  'feu 
W"(u)  =  7 


Setting 


we  obtain 


W'(  uj)  =  D  +  70;  =  t/x 


t  —  Dx 


LUq  = 


7X 


(770) 

(771) 


(772) 

(773) 

(774) 

(775) 
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Now, 


and  also 


In  addition 


where 


K(u)  =  Du  +  ^u2  /  2 

1 


2^x 


( t  —  Dx)(t  +  Dx) 


ust  —  K{us)x  = 


(t  —  Dx)* 
2'yx 


us  —  u0  =  — (t  —  Dx  —  Uo'yx) 
'yx 

=  -(t-(t)x) 

7 x 


(t)x  =  (D  +  Uo'y)  x 


(776) 

(777) 

(778) 

(779) 

(780) 

(781) 


The  asymptotic  solution  is 


ua(t,x)  ~  F(u>s,  0)  Y  “e  l7r/4exp 


(t  —  Dx)* 
i - 

2^x 


(782) 


But 


(a/7r)1//4 

F(u,  0)  = - — —  exp 


Vv 

(a/n  )1//4 
\Ja  —  i/3 


exp 


(t  —  Dx  —  Uo'yx)* 

2r/'y2x2 

aft  —  Dx  —  u0yx)2  Pit  —  Dx  —  u^x)2 

i- 


2(a 2  +  P2)  72x2  2(a2  +  P2)^2x2 


(783) 

(784) 


and  therefore 


/.  v  (a/7r)1/4  [T  _i7r/4 


a/o  —  ip  V  7^ 


a(£  —  Dx  —  Uo'yx)2  .pit  —  Dx  —  Uq^x)2  +  .  (t  —  Dx)* 


2(a2  +  P2)  *y2x2  2  (a2  +  P2)^2x2 


\ua(t,x)\2  = 


1  (a/ 7T)1/2 


7^  \/ a2  +  /52 


exp 


2  i 


Oi{t  -  (t  )x) 


(a2  +  P2)  72x2 


2'yx 

(785) 

(786) 


At-(t)x)2 

2rp2x2 

At-(t)x)2 

2'p2x2 


1 

7a; 


(f  —  Ha;)2 
2jx 

~ — 2^^  “ — ^  wo(*  —  ( t  )x)  + 

+Mt~(t)x)  +  ^2'yx 


(787) 

(788) 

(789) 
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where  we  have  used 


(790) 

(791) 


t  —  Dx  —  t  —  (t)x  +  ca0  'yx 
( t  -  Dx)2  =  (t-  (t  )x)2  +  2{t  -  ( t  )x)u0'yx  +  affix2 


11  Future  Research 

There  are  a  number  of  directions  that  the  above  idea  of  pulse  propagation  should  be  extended 
to. 

1.  The  case  with  damping  should  be  investigated  and  the  formalism  developed  here 
should  be  extended  to  handle  that  case. 

2.  Can  these  methods  presented  be  generalized  to  equations  with  non-constant  coeffi¬ 
cients? 

3.  The  Gabor  procedure  for  a  pulse  needs  investigation.  That  is  a  proper  analysis  of  the 
concept  of  an  analytic  signal  for  a  pulse  has  not  been  done. 

4.  Instantaneous  frequency  and  local  spatial  frequency  have  been  obtained  for  the  as¬ 
ymptotic  solution.  One  should  try  to  improve  on  this. 

5.  Exact  calculation  of  moments.  We  have  shown  that  the  moments  can  be  calculated 
exactly  and  easily.  Can  these  moments  be  used  to  construct  a  better  approximation  than  the 
classical  asymptotic  approximation?  That  is  construct  an  approximation  to  \u(x,  t)\2  and  in 
particular 


Find  an  approximation  to:  \u(x,t)\2  Given:  (792) 

u(x, 0)  and  (793) 

(xn)t  n  —  1,  N  (794) 

There  are  many  methods  to  construct  densities  from  a  given  set  of  moments  and  it  would  be 
interesting  to  apply  these  methods.  Note  that  in  this  formulation  one  would  only  approxi¬ 
mate  \u(x,  t)\2  ,  that  is  the  magnitude.  But  if  we  add  to  this  the  time  moments,  then  it  may 
be  possible  to  also  get  the  phase  to  a  better  approximation. 

6.  Can  the  methods  and  models  developed  here  by  applied  to  nonlinear  wave  equations? 

7.  The  accuracy  of  the  Wigner  approximation  scheme  needs  further  investigation. 

8.  Can  one  obtain  equations  of  motion  for  the  amplitude  and  phase  separately?  We  write 

u(x,t)  =  Rixrfe**^  (795) 
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We  have  been  able  to  show  that 

-7^-  =  —  Im  u*  (x,  t)W ( KI)u(x ,  t )  (796) 

-7^-  =  — Re  «*(;£,  t)W(tC)u(x,  t)  (797) 

Can  these  equations  be  solved  directly  for  phase  and  amplitude? 

9.  The  Wigner  distribution  approach  should  be  generalized  to  other  distributions  \cite{cohen66, rev, bool 
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A  Appendix:  Notation  for  Part  A 


k 

UJ 

u(x,t)  =  A(x,  t)el^x,t^ 

A(x,  t),  (p(x,  t) 
u  =  W(k ) 
v(k)  =  W'(k) 

u(x,t)  =  -^=  f  S(k,  0)  elkx-'lW(k)t  dk 
S(k,  0)  =  -j=  f  u(x,  0)  e~lkx  dx 
S{k,t )  =  S(M) 

=  -A=  f  u(x,  t)  e~lkx  dx 

S(k,t )  =B(k,t)e^k^ 

B(k,t),  ip(k,  t) 

x  =  iK 

( xn)t  =  f  xn  \u(x,t)\2  dx 
V  =  f  v(k)  \S(k,0)\2dk 

Covxt,|j 

P  ^ 

tc 

tM 

Ui(x,t)  =  --§-t<p{x,t) 
ki(x,t )  =  J^OM) 


Spatial  Frequency  (wave  number) 

Frequency 

Pnlse  at  position  x  and  time  t 
Amplitude  and  phase  of  pnlse 
Dispersion  relation 
Group  velocity 
General  solution  for  a  pnlse 
Spatial  spectrum  at  time  zero 
Spatial  spectrum  at  time  t 
Spatial  spectrum  at  time  t 

Spatial  spectrum  in  terms  of  amplitude  and  phase 

Amplitude  and  phase  of  spatial  spectrum 

position  operator  in  k  space 

Spatial  moments  of  a  pnlse  at  time  t 

Average  group  velocity 

Standard  deviation  of  group  velocity 

Covariance  of  position  and  group  velocity  at  time  t 

Correlation  coefficient 

Time  of  contraction 

Time  at  minimum  contraction 

Instantaneous  frequency 

Local  spatial  frequency 


79 


B  Appendix:  Notation  for  Part  B 


k 

C 0 

u(x,t)  = 


A(x,t), 

k  =  K( 

uj) 

t(uj)  = 

K\u) 

u(x, t ) 

=  e<wt~ 

iK(i 

F(k,  0) 

=  7^  fu(x’°)  e~iUt 

dt 

F(u,  x) 

=  F(k,  0)  e~iK^x 

F(u,  x) 

=  vb  fu(x >*)  e~iWt 

dt 

F(u,  x) 

=  B(  u,x)e>^u^ 

B(u,  x) 

T  =  i-g- 

OUJ 

(tn)x  = 

=  f  tn  \u(x,  t)  2  dt 

T  =  f  r(u)  \F(k,0)\2  du 

<?T 

Covj-j-j^ 

Pt.T  \x 

Xc 

Xm 

Ui(x,t) 

= 

ki(x,t) 

=  -mvM 

Spatial  Frequency  (wave  number) 
Frequency 

Pulse  at  position  x  and  time  t 
Amplitude  and  phase  of  pulse 
Dispersion  relation 
Transit  time 

General  solution  for  a  pulse 

Spectrum  at  time  zero 

Sspectrum  at  position  x 

Spectrum  at  position  x 

Spectrum  in  terms  of  amplitude  and  phase 

Amplitude  and  phase  of  spectrum 

Time  operator  in  u  space 

Time  Moments  of  a  pulse  at  x 

Mean  transit  time 

Standard  deviation  of  transit  time 

Covariance  of  time  and  transit  time  at  x 

Correlation  coefficient 

Positions  where  pulse  is  contracted 

Position  of  minimum  contraction 

Instantaneous  (time)  frequency 

Local  spatial  frequency 
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C  Appendix:  Moments 


In  this  appendix  we  give  the  derivations  of  the  moments  of  a  pulse.  As  pointed  out  in  the 
text  the  spatial  moments  of  a  pulse  as  a  function  of  time  can  be  calculated  directly  and 
easily  from  the  initial  spectrum.  In  particular  we  have  that  in  general: 


( xn)t=  /  xn  \u(x,  t)\2  dx 


=  /  S*(k,t)  Xn  S(k,t)  dk 


where  X  is  the  position  operator  in  the  k  representation 


X 


d_ 

dk 


(798) 

(799) 


(800) 


And  again  as  we  pointed  out,  what  makes  the  calculation  easy  in  the  spectral  Fourier 
domain  is  that  the  time  dependent  spectrum  is 


S(k,t)  =  S(k,  0)e~lW{k)t 


(801) 


First  we  list  here  different  expressions  for  the  moments  that  are  useful 


(x)t=  1  x  \u(x,  t)  2  dx 

(802) 

=  1  S*{k,t)XS(k,t)dk 

(803) 

x2)t=  j  x2  \u(x,  t)  2  dx 

(804) 

=  1  S*(k,t)X2S(k,t)dk 

(805) 

=  j  \XS(k,t)\2dk 

(806) 

ax\t  =  J (x  -  (x  )t)2  \u(x,  t)\2  dx 

(807) 

=  (x2)t-((x)t)2 

(808) 

=  1  S*(k,t)(X  -  (x)t)2S(k,t)dk 

(809) 

=  1  \(X-(x)t)S(k,t)\2dk 

(810) 

First  Moment.  We  have 


(x)t 


x  | u(x,  t)\2  dx 


S*(k,t)XS(k,t )  dk 


(811) 
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but 


and  therefore 


XS(k,t)  =  i-^-S(k,t)  = 
ok 


i-^-S(k,  0)  +  tvS(k,  0)  )  e~iW{k)t 
ok 


{x)t=  I  S*(k,0)XS(k,0)dk  +  t  J  v\S(k,0)\2dk 
—  (x)0  +  t  j  v\S(k,  0)|2  dk 
=  (x)0  +  t(v)0 


(812) 

(813) 

(814) 

(815) 


Second  Moment.  For  the  second  moment  we  do  it  two  different  ways.  First  by  way  of  Eq. 
(806)  and  then  directly  by  way  of  Eq.  (805).  We  calculate 


\XS(k,t)\2  =  I  [i-^S(k,0)+tvS(k,0)  )  e-iW^ 

iSrS(k,  0)  +  tvS(k ,  0) 

Ok 

=  [-i-^S*(k,0)  +  tvS*(k,0)^j  (i^S{k,Q)+tvS{k,Q) 

d 


(816) 

(817) 

(818) 


d 


=  | XS(k,  0)|2  +  t2  v2 1 S(k,  0)|2  +  f  t[vS*(k,  0 )—S(k,  0)  -  vS(k,  0 )—S*k(k,  0)] 


dk 


dk 


and  therefore 


( x2)t=  /  \XS(k,  0)|2  dk  +  t2  /  v2  \S(k,0)\2  dk 


=  it 

=  it 


K*(M)J)s(M)  -  vS(k,o)^s;(k,o) 

r\  r\ 

vS*(k,  Q)—S(k,  0 )+S*{k,  0)—vS(k,  0) 


dk 


where 


=  {x2)o  +  t  J  S*(k,  0)[u,  X]+S(k,  0)  dk  +  t2  ( v2  ) 
=  { x2  )0  +  t  ( [v,  X]+  )0  +  t2{v2) 


r  ,  d  .  d 
[v,X]+=,v-+l-kv 

=  vX  +  Xv 


(819) 


(820) 

(821) 

(822) 

(823) 

(824) 

(825) 
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We  now  do  it  the  second  way,  that  is  by  way  of  Eq.  (805), 

<92 

X2S(k,t)=i2—S(k,t) 

=  t2-^-S(k,  0)  e~iW^ 
ok 2 

-kLS(k,  0)  +  it  S(k ,  0)  +  2v^)  +  vVS(k,  0) 


D—iW(k)t 


and  therefore  we  obtain 

(x2)t  =  (x2)0  +  it  J  v  ^S*(k,0)^S(k,0)k  -  S(k,0 dk  +  t2(v2)0 

=  (x2)0  +  it  J  S*(k,  0)  0)  dk  +  t2(v2)0 

=  {x2)0  +  t  { [v,X]+)o  +  t2  (v2) 

as  before. 

Standard  Deviation  or  Spread.  The  standard  deviation  is  therefore 

<*l\ t  =  {x2)t-  {x)2  =  (x2)0  +  t2  (v2)  +t([v,X]+)0-[(x)  +t(v)  }2 


or 


where 


<%\t  =  ax\o  +  2t  Covxv  +  t2  a2 


Cov(vx)  =  -(tW  +  Xv)q  —  (v)o{x)o 


Global  k  moments.  We  have 

(k)t  =  f  k\S(k,t)\2dk 


=  J  u*(x,  t)tCu(x,  t)  dx 
(  k2)t=  j  k2  \S(k,t)\2  dk 

=  /«•(*,  & 

=  j  \JCu(x,  t)\2  dx 

<?k\ t=  j{k-(k)t)2\S(k1t)\2dk 
=  (k2)t~((k)t)2 
=  j  u*(x,t)(lC  —  (k)t)2u(x,t)  dx 

=  (  |(/C  -  (  k)t)2)u(x,t)\2  dx 


(826) 

(827) 

(828) 

(829) 

(830) 

(831) 

(832) 

(833) 

(834) 

(835) 

(836) 

(837) 

(838) 

(839) 

(840) 

(841) 

(842) 

(843) 
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In  the  case  of  wave  propagation  these  moments  are  constant  in  time. 

Covariance.  To  calculate  the  time  dependence  of  the  covariance  between  position  and  wave 
number  we  first  note  some  general  properties  of  the  anticommutator  operator  when  the 
variables  are  x,  k.  Using  the  commutator  relation 

[X,K\=i  (844) 

we  have  that 

[X,  JC]+  =  2ICX  +  i  —  2XIC  —  i  (845) 

Also  for  two  arbitrary  functions  a  and  b 

[X,  K]+ab  =  a[X,  lC\+b  +  b[X,  K]+a  —  ab  =  2 kXab  +  iab  (846) 

Now  applying  this  to  S(k,  0)  we  have  that 


(847) 


[X,  JC\+S(k,  0)  e~iW{k)t  =  2  ik(Sk(k,  0  )-ivt  S(k,  0))  e~iW{k)t  +  S(k,  0)  e~iW{k)t  (848) 

=  e~iW{k)t[X,  JC]+S(k,  0)  +  2 kv  tS(k,  0)  e~iW{k)t  (849) 


Hence 


{xk)t  =  {-[X,K]+)=  /  -S*(k,0)[X,K]S(k,0)dk  +  t  /  kv\S(k,0))\2dk 


or 

{xk)t  —  (xk)0  +  t(xv)o 


(850) 

(851) 


Thus  ( xk)t  increases  linearly.  Using  the  fact  that  all  wave  number  averages  are  independent 
of  time  we  have 


Co vxk(t)  =  (xk)t-  (x)t(k)t  (852) 

=  (xk)t-  (x)t(k) o  (853) 

=  (xk)0  +  t  (xv)0  -  (k  )0((  x  )0  +  t{  v ))  (854) 

or 

Co vxk(t)  =  Covxk(0)  +  t  Covxv(0)  (855) 

The  Covariance 

We  now  consider  the  covariance  between  x  and  k: 

Covxk\t  =  (xk)t  -  (x  )0(fc)o  (856) 
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We  have 

(xk)t  =  l(  kX  +  Xk  )t  =  [  k^^hsfct)]2  dk 

J  ok 

(857) 

and 

(x  k)t  =  (x  k)  o  +  (k  v  )t 

(858) 

Inserting  this  into  Eq.  (855)  gives 

Govxk\t  Cov^jo  “I-  t  Cov kv 

(859) 

where 

Covkv  =  (kv)  -  (k)(v) 

(860) 

We  note  that  Covkv  is  independent  of  time. 

Also,  one 

can  obtain  that 

+  Xv)t  =  \(vX  +  Xv)o  +  {v2(k)  )t 

(861) 

and  hence 

CoV 1 1  Covxfc|0  H“  &vt 

(862) 

_  Cov AK\t  _  Co vxk(t)  =  Cov^O)  +  alt 

pxk\t  j 

ax\tak\t  av  j (j2|0  +  2 1  Covxv  +  t2al 

(863) 

As  before  we 

have  that 

Pxv\t  *  1  ^ 

(864) 

D  Appendix:  Marginals  of  Space-Time  Distributions 

We  give  the  marginals  of  the  four  dimensional  distribution  given  by  Eq.  (383)  [7].  First,  for 
convenience  we  repeat  the  fundamental  definition 


u(x,  t ) 
S(k,t) 

F(u,  x ) 

G(k,u) 


A{x,t)eiipixA 


(865) 

(866) 

(867) 

(868) 
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and  also  repeat  the  four  different  ways  one  can  write  the  four  dimensional  distribution 


2 

1 

-t 

2 

¥ 

1 

'  2 


1 

2 ' 


D~iTxk—iru} 


dr  dT,r 


1 

-( 

2 


1 

2' 


1 

-( 

2 


1 

-  -( 
2 

¥■ 

1 

-  -T 

2 


1 

-t 

2 


1 

2 


1 

-t 

2 


There  are  4  three  dimensional  marginals 


(869) 

(870) 

(871) 

(872) 


P(x,k,t)  =  J  C(x,  k,t,co)  du 

1  r  i  x 

=  —  /  u*(x  -  -rx,  t )  m(x  +  -r*,  t)  e~JTxk  drx 

=h  !  s*(k+l6^t)s(k~ld^t)e~j9xXdd- 


P(x,k,co)  =  /  C(x,k,t,u)dt 


P(x,  t ,  oo) 


P(k,  t ,  oo) 


1  G*(k  +  ^9x,co)G(k-^9x,u)  <  dO., 
j  F*(x  -  ^rx,u)F(x  +  ^Tx,u)e~JTxk  drx 
J  C(x,  k ,  t,  oo )  dk 

J  u*(x,t-^r)u{x,t+^T)e~JTUJ  dr 

f  F*(x,oo  +l-6)F{x,u -\Q)e-*»  d9 


C(x,  k,  t,  o o)  dx 


1 

27T 

1 

27T 


G*  (k,  to  +  U)  G(k,  oo  —  U)  e-j6t  d,9 
S*(k,  t  -  ir)  S(k,  t  -  ir)  e~jTU}  d9 


(873) 

(874) 

(875) 

(876) 

(877) 

(878) 

(879) 

(880) 
(881) 
(882) 

(883) 

(884) 


There  are  six  two  dimensional  marginals  given  MUST  BE  CHECKED 


PM=  /  W{x,  k,  t,  u)  dx  dk 


—  —  if  u*(x,t  —  -t)  u(x,t  +  -t)  e  lrud,Tdx 
27 t  J  J  2  2 

P(x,k)—  j  W(x,k,t,u>)  dtdcu 

=  J  J  U*(x-  ^tx,  t )  u(x  +  ^tx,  t )  e~tTxk  drxdt 

P(k,  t)  =  j  W (x,  k,  t,  u)  du)dk=  \S(k,t)  \ 2 
P(k,  u)=  j  W(x,k,t,u))dtdx—  \G(k,  cu)  |2 
P(t,  x)  —  j  W(x,k,t,u>)  du dk=  \u(x,t)  |2 
P(x,  u)=  f  W(x,  k,t,u>)  dtdk—  \F(u,x)\2 


There  are  4  one  dimensional  marginals 


P(t)=  j  W(x,  k,t,u)  dx  dkdu  — 
P{x)—  j  W(x,k,t,u)  dtdudk  = 
P(k)  —  j  W(x,  k,t,u>)  dcu  dkdt  = 
P{ u)=  (  W(x,k,t,u>)  dtdkdt  —  I 


I  |  u(x,  t )  1 1dx 

r 

|  u(x,  t )  1 2dt 
(S^/c,  t)  | 2  dt 
| F(u,  x)  |2  dx 


(885) 

(886) 

(887) 

(888) 

(889) 

(890) 

(891) 

(892) 

(893) 

(894) 

(895) 

(896) 
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Frequency  moments.  The  first  conditional  moment  of  frequency  is 


{v)x,k,t  =  /  u>C(x,k,t,u>)  du 


(897) 


1  1 

V - Tx,  t - ' 

2  2 

1  1  \ 
t--T)u\ 


1  1 

^  +  +  - 

1  1 


—irxk 


-ir^-iru  drdTxdw  (ggg) 

i  drx  (899) 


1  1 
27 r  2  i 

1  1 
27 r  2i 


Ft 


d 


U*(X  -  -Tx,t)u(x  +  - Tx,t )  ~U*(X-  - Tx,t )  — «(£  +  ifx.t) 


„  — irxk 


d.Tx 


(900) 


1  x  d  ..  1  ,  1  sd  1 

A^x+2Tx,t'dtA^X^  2Tx,t'  ~  A^X  ~  2Tx,t'FA^X+  2Tx,t' 


1  \  A  f  1  J  f  ,  1  \  ,  1 

-lA[X  +  -Tx,  t)A(x  -  ~TX,  t)—  {  (p{x  +  —Tx,  t)  +  v(x  -  -Tx,  t) 

ei<t>{x+±Tx,t)-i<j>(x-±Tx,t)-iTxk 


Also, 


Using  Eq.  (897)  we  have  that 


{w)x,t  =  /  n >C(x,  k,  t,  lu)  duidk 


(901) 


(902) 


J  uC(x,  k,  t,  u)  dudk 

1 


"  d_ 

2m  J  d r 

i  r 


d 

'a? 


1  1  N 

-rx,  t - r)n 

2  2  ’ 

1  1 

7~T  1  t  ^ 
2  2 


1 

2 


—irxk 


1 

2 


1 

2 ' 


*  J  dr 

which  evaluates  to 


d  .  .  1  ,  .  1 

U  (X,t-  -T)u(x,t+  -T)\r=0dTx 


9  <9 

=  A2(x,t)—if(x,t) 


(903) 

r=o 

(904) 

drx 

(905) 

(906) 

(907) 

Conditional  Position.  The  average  position  is  similarly  given 


( x)k,t,u>=  /  xC(x,  k,t,u)  dx 


(908) 


■  ii  I K-  r>  r>  !•-“* 


1  1  f  l_^d 'B{k,t  +  \r)  1  dB(k,  t  —  \r) 


.  B(k,t - r) 

27r  2i ./  v  ’  2  J  dk 


S(M+W  dk 


(909) 

(910) 


1  1  Q  C  ~\  1 

+  jB(k,  t  +  -r)B{k,  t  -  -t)—  \  ip{k,  t  +  -r)+  ip(k,  t--r)  \ 

2  2  (jrC  I  2  2 


To  obtain  the  average  position  for  given  wave  number  and  time  we  have 

{x)k,t  —  J  xC(x,k,t,cu)  dxdu 

=  B2(k,t)  J^(M) 


We  now  further  average  to  obtain 


d 


d 


(x)t=  J  B  ( k,t )  t)  dk  —  { t ) ) 

Covariance.  We  now  calculate  the  covariance.  Consider  ( ku))x>t 

(ku)xj=  j  ku>C(x,k,t,u>)  dkdu 

=  j  k{u>)Xtkjt  dk 

a2(  .s& <p(x,t)  d <p(x,t) 

^  ’  >  dx  dt 

The  local  covariance  of  frequency  and  spatial  frequency  is  therefore 

'  d(p(x,t)  d<p(x,t)  d<p(x,t)  d(p(x,t) 


CovXtt(ku>)  =  A2(x,t ) 
and  the  global  covariance  is  hence 


dt 


dx 


dt 


dx 


(911) 

(912) 

(913) 

(914) 


(915) 

(916) 

(917) 

(918) 


Co  v(kuj)  — 


A2,x  (  d<p(x,t)  dip(x,t)  _  dip(x,t)  _  d(p(x,ty  (  ^  ^ 


dt 


dx 


dt 


dx 


(919) 


As  discussed  previously,  we  can  think  of  and  as  the  frequency  and  spatial 

frequency  in  the  x,  t  representation.  Using  this  idea  for  example  we  can  immediately  write 


{k)x,t  =  A2(x,t) 


d(p(x,  t) 
dx 


In  the  spectral  domain  we  have 


{xt)k,  oj  =  B2(k,u ) 


dip(k,u)  &ip(k,u) ) 


d  k 


du 


and  again  one  can  think  of  x,t  as  ■ 


(920) 


(921) 


E  Appendix:  Asymptotics  for  pulse  propagation 


In  the  text  we  have  given  formulas  for  the  instantaneous  frequency  where  we  have  based  our 
approach  on  the  classical  asymptotic  approximation.  We  repeat  here  the  classical  asymptotic 
result.  The  pulse  is  given  by  (exactly) 


u(x,  t ) 


(922) 


where  S(k,  0)  is  the  initial  spatial  spectrum  and  it  is  calculated  from  the  initial  pulse.  The 
asymptotic  solution  as  standardly  given  is 


ua(x,t)  ~  S(ks,  0)< 


tW"(ks 


0iksx— iW(ks)t— i7rsgn  W"  /  4 


(923) 


where  k*  is  obtained  from 


x 


W\ks)  =  - 


In  deriving  the  formulas  for  instantaneous  frequency  we  wrote 


ua(x,t)  B(k,  0)e 


iip(k,  0) 


tW"(ks 


0iksx—iW (k3)t—iiTsgn  W" / 4 


(924) 


(925) 


where  B(k,  0)  and  ijj(k,0)  are  the  initial  amplitude  and  phase  of  the  spectrum.  The  asymp¬ 
totic  amplitude  and  phase  of  the  pulse  are  therefore 


U*,t)  =  \s(k',  o)l^^y 

<Pa(x,t)  =  Ip{ks,  0)  +  ksx  -  W(ks)t  -  7rsgn  W" / 4 


(926) 

(927) 
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and  we  obtained  onr  formulas  for  instantaneous  frequency  by  differentiating  asymptotic 
phase  appropriately.  We  repeat  those  results  here 


k  (x  t)  =  °(x  t )  =  1  gfhh 

gxif(x,  )  tw„^  dk* 


+  ks 


(928) 

(929) 


We  now  show  that  a  different  approach  is  more  appropriate  and  gives  considerable  better 
answer.  This  has  recently  been  done  with  P.  Loughlin  [14].  Our  approach  consists  of 
two  steps.  First  we  derive  a  more  suitable  asymptotic  approximation  then  the  classical 
one.  Secondly,  using  this  new  asymptotic  approximation  we  calculate  the  Instantaneous 
frequency. 

To  derive  the  new  aproximation  we  take  the  initial  spectral  phase  into  account  in  deriving 
the  asymptotic  formula.  We  write 


u(x, t ) 


B(k,0) 


i kx—iW(k)t+iip(k,0 ) 


(930) 


In  the  stationary  phase  method  one  assumes  that  js  rapidly  oscillating  and 

the  amplitude  is  relatively  slowly  varying.  Hence,  for  most  regions  of  integration  there  will 
be  cancellations.  However  if  there  is  a  region  where  the  phase  is  not  oscillating  then  that 
region  will  contribute  the  greatest  part  to  the  integral.  That  region  is  where  the  derivative 
of  the  phase  is  zero 

d 

-[kx-W(k)t  +  ip(k,  0)]  =  0  (931) 

OK 


This  gives 


x-W\k)t  +  ^\k,  0)  =  0 


(932) 


or  equivalently, 

w\k) = 

One  solves  this  equation  for  k  and  writes  the  solution  as  ks.  Notice  now  that  ks  is  no  longer 
a  function  of  x/t  as  would  be  the  case  if  we  did  not  take  into  account  the  phase  of  the  initial 
spectrum.  Since  we  assume  that  the  amplitude  is  slowly  varying  we  have 

ua(x,t) - ^B(ks,  0)  [  e*x-iW(kyt+Wk,o)  dk  (934) 

V  2vr  J 

The  phase  is  now  expanded  in  a  series  about  the  stationary  point  ks  and  one  keeps  terms 
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only  up  to  the  quadratic  one.  In  particular 


kx  —  W ( k)t  +  ip(k,  0)  ~  ksx  —  W ( ks)t  +  ip(ks,  0)  +  [x  —  W'{ks)t  +  ip'(ks,  0)]  (k  —  ks )  (935) 

+  \[-W“{ka)t  +  i>"(ks,  0 )](k  -  ks )2  (936) 

=  ksx  -  W (ks)t  +  ip(ks,  0)  +  ^[~W"(ks)t  +  ip"(ks,0)](k  -  ks) 2 

(937) 


and  therefore 


ua(x,t) 


1 


B(ks,  0)  /  exp  i[ksx  -  W(ks)t  +  ip(ks,  0)  +  -[—W"(ks)t  +  ip"(ks,  0)](/c  -  ks )2]  dk 


V7^ 


S(ks,  0)e 


(938) 

expi[-{—  W"{ks)t  +  ip"(ks,  0)}(A;  —  ks )2  ]  d,k  (939) 


Using, 


e^2/2cft  = 


27T  /  27T  • 


jTT/ia/4 


-ja  V  a 


where  /ia  is  the  sign  of  a 


Ha  =  sgn(a) 


(940) 


(941) 


we  have 
ua(x,t) 


V7^ 

S(ka,  0) 


S(ks,  0)e 


i[fcs:r— w(fcs)t] 


exp  i[-{-fU"(A;s)t  +  ^"(^s,  0 )}(k  -  ks  )2  ]  dk  (942) 


tW"(ks)  -  ip" (ks,  0) 


0iksx— iW(ks)t—  msgn  [W"t— ■0,/  (A:s  ,0)]  / 4 


(943) 


Note  the  following.  From  a  functional  point  of  view  there  appears  to  be  only  a  difference 
in  the  amplitude  and  not  the  phase  when  one  compares  Eq.(943  with  Eq.  (925).  However 
it  must  be  emphasized  that  there  is  a  difference  as  to  the  value  of  the  stationary  points. 
Namely,  for  this  case  we  must  solve  Eq.  (  933)  rather  than  Eq.  (924)  for  ks. 

Instantaneous  frequency.  We  now  obtain  the  instantaneous  frequency.  The  phase  is 


(pa(x,t)  =  ip(ks,  0)  +  ksx  -  W(ks)t  -  7Tsgn  [W"t  -  if"(ks,  0)]/4  (944) 
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Differentiating  the  phase,  (pa(x,t),  we  have 


.  .  d 

Vi(X,t)  =  <pa(x,t ) 


dip  dks  dks  dW (ks)  dks 

dk.  dt  X  dt 


dk,  dt 


+  W(ks 


dip  |  x  tdW(ks 

dks  dkg 


w+w^ 


dv  ( x  +  i/(/.\  0) 


dkg 

dt 


+  W(ks 


But  since 


W\k)  = 

the  first  term  is  identically  zero  and  hence 


x  +  ip'(k,  0) 


Ui(x,t)  =  W(ks) 


(945) 

(946) 

(947) 

(948) 

(949) 

(950) 


Also, 


d 

ki(x,t )  =  —tpa(x,t) 
dx 

dip  dkg 


_ +  x^  +  k  dW{kg)dks 

dkg  dx  dx  s  dkg  dx 


<kp_  _  dW(k 
dkg  +  dkg 


dkg 


dib 

dk.+x~t 


dx 

x  +  ip'(k,  0) 
t 


+  ks 


dkg 

~a - ^  ^  s 

dx 


and  therefore 


h(x,t)  =  kg 


(951) 

(952) 

(953) 

(954) 

(955) 


Example.  We  now  take  an  example  to  illustrate  the  differences.  We  shall  do  it  three 
different  ways,  exact,  using  the  classical  asymptotic  approach,  and  the  modified  asymptotic 
approach  discussed  above.  At  t  =  0  we  take 


u(x,  0)  =  eifSx2/2 


and  the  initial  spectrum  is 


S{k,Q) 


For  the  dispersion  relation  we  take 


W(k)  =  ,yk2 /2 


(956) 


(957) 


(958) 
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Exact.  The  time  dependent  spectrum  is 


/  2 

k2 

rw 

-i^-iw(k)i 

giving 


u(x,  t )  = 


S(k,  t)e 


ikx 


i  i  r 

yfhtyJ-iP  J 

1 

.  =  exp 

a/1  +  7  f3t 


exp 


.  k 2  i'ykH 

-'w~—+  x. 


dk 


,/3  x2 


The  exact  phase  is  therefore 


ip{x,t)  = 


2  1  +  7 /?*. 
(3  x2 


2  1  +  7  Pt 


and  the  exact  instantaneous  frequency  is 


and  the  exact  spatial  instantaneous  frequency  is 

Standard  asymptotic  method.  The  stationary  points  are  obtained  from 

x 


and  hence, 


giving 


W\k)  =  7  k  =  - 


ks  =  — 

7t 


W(ks)  = 


x 


2  7t2 
x2 


ksx  -  W{ks)t  =  — 


Therefore, 


ua(x,t)  ~  S(ks,  0/ 


tW"(ks 


0iksx—iW (ks)t-iirsgn  K" /4 


7  Pt 


exp 


.  x  .  x 
2/?72t2  27 1 


(959) 

(960) 

(961) 

(962) 

(963) 

(964) 

(965) 

(966) 

(967) 

(968) 

(969) 

(970) 

(971) 
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The  phase  is 


r2  2 
^(X'*)=  2fHV  +  2it 

and  differentiating  we  obtain 

r2  2 

(972) 

(973) 

As  a  check  we 

also  calculate  the  instantaneous  frequencies  by 

(974) 

But 

dfi  k  x 

dks  fi  fi'yt 

(975) 

and  therefore 

2 

<"y>  nr>  ^ 

/  \  tX/  tv  tv 

k’i(X'()=  2-yt* 

9  9 

(976) 

X  X 

/?72f3  '  27f2 

(977) 

Also, 

,  ,  .  1  dfi 

k‘M=  tw"(k,)dk,+k- 

(978) 

X  X 

fi^H2  7 1 

(979) 

As  can  be  seen  these  are  approximations  to  the  exact  answers. 


Modified  Asymptotic  method.  We  have  to  solve 


W'{k)  =  k 7 


x  +  ip'(k,  0) 
t 


That  is, 


and  therefore 


x  -  ks/fi 
hi  =  — — 


(980) 

(981) 


_  fix 

s  1  +  7  fit 

(Notice  that  this  happens  to  be  the  exact  instantaneous  spatial  frequency,  but  we  do  not 
take  that  into  account.) 
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First,  we  calculate  W(ks)  and  ksx  —  W(ks)t.  In  particular 

W(k9)  =  7^/2 

fix 


7 

2  VI +  7 fit 

7  fi2x2 

2(1  +  7  fitf 


and  also 


ksx  —  W  ( ks)t  = 


fix 


7  fi2x2 


-x 


1  +  7  fit  2(1  +  7  fit)2 

fix2  f,  7  fit 

1  +  7  fit 
fix2 


t 


1  - 


(1  +  7  fit)2 


2(1  +  7  fit) 
{l+lfit/2) 


Therefore 

ua(x,t)  ~  S(ks,0)i 

■  exp 


tw"(ks)-r(K  o) 

.  k2 


0iksx—iW (k„)t—iirsgn  [W"t—  ■0,/ (fcs ,0)] / 4  _ 


—  I 


2  fi. 


fi  V  t'y  +  1  /fi 


exp 


tW"(ks )  -  'fi'fiks,  0) 
fix2 


vikax—iW (k3)t—iirsgn  [W"t—il>"(ka, 0)]/4 


fi  V  *7  +  1//5 


exp 


'(1  +  7/^)2 

.  /la:2 

'(1  +  7  fit)2 

fix2 

fi  V  <7  +  1//3  L’  (!  +  7W 

fix2 


(1  +7^/2)  -  + 


fix 


2/1  V 1  +  7/11 


(1  +7/3f/2)  -1- 


fix 


■  exp 


2/1  V 1  +  7  fit 

(1  +  7/11/2-1/2) 


1  +  I7/I 


exp 


2(1  +  7  fit) 


(983) 

(984) 

(985) 

(986) 

(987) 

(988) 

(989) 

(990) 

(991) 

(992) 

(993) 

(994) 


But  this  is  the  exact  answer  and  hence  will  give  the  exact  instantaneous  frequencies.  To 
verify  we  use 

7  fi2x2 


and 


Vi(x,t)  =  W(ks )  = 
ki{x,  t)  =  ks  = 


2(l  +  7/ll)2 
fix 


(995) 


(996) 


1  +  7  fit 

and  indeed  these  are  exact.  We  therefore  see  that  the  approximation  we  have  presented 
here  is  more  accurate  than  the  standard  asymptotic  one. 
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Case  B 


We  write  the  corresponding  equations  for  case  B.  The  general  solution  is 


u(x,t )  = 


V2k  . 


F(u,x)  = 


F(u,x)  e  du 


u(x,t)e  luJt  dt 


with 


F(u,x)  =  F(u,  0)  e~iK^)x 

and  where  F(u,  0)  is  the  spectrum  evaluated  at  x  —  0, 


F(u,0)  = 


The  standard  asymptotic  solution  is 


u( 0,  t)  e  lujt  dt 


ua(x,t)  ~  F(u,  0)< 


0iudt— iK(u)x— 27rsgn  K"  /  4 1 


where  one  obtains  us  from 


The  amplitude  and  phase  are 


xK"(u ) 
K'iojg)  =  t/x 


\UJ=(jJs 


Mm)  I  =  I^K°)K,xX„(w) 

<pa(x,  t)  =  0)  —  K(u)x  +  ut  —  7TSgn  K " /A 

and  the  instantaneous  frequencies  are 

.  .  1  di/j 

=  77777  ^—+US 


h(x,t)  = 


xK"{u)  du 
t  difj 


+  K  (c os 


x2K"(us)  dus 

We  no  obtain  the  new  formulas.  Using  the  same  approach  as  before  we  now 

u(x,  t )  =  -L  [  B(u,  0)  e^t-iK(u)x+i^u,o)  du 
v27 X  J 

and  the  region  where  the  derivative  of  the  phase  is  zero  is  obtained  from 

d 


du 


[ut  —  K{ u)x  +  0)]  =  0 


(997) 

(998) 

(999) 

(1000) 

(1001) 

(1002) 

(1003) 

(1004) 

(1005) 

(1006) 
have  write 

(1007) 

(1008) 


97 


which  gives 


t  —  K'(uj)x  +  ip'(u),  0)  =  0 


(1009) 


K'(u)  = 


t  +  ip\u,  0) 


Taking  the  amplitude  out  of  the  integral  we  have 


ua(t,x)  ~  -=£(u;s,  0)  /  e^-iK(u,)x+i^,o)  du 
V  2vr  J 


(1010) 


(1011) 


Expanding  the  phase  about  the  stationary  point  us  and  keeping  terms  only  up  to  the 
quadratic  one  we  have 

tut  —  K(u)x  +  ip(u,  0)  ~  u jst  —  K(ous)x+ip(ous,  0)  +  -[— K"  (us)x  +  ip"  (us,  0)](o;  —  us)2  (1012) 
and  therefore 

If  1 

ua(t,x)  ~  -j=B{us,  0)  J  expi[ujst- K(ujs)x  +  ip(ujs,0)  + -[-K"(lus)x  +  iP"(lus,0)](uj -ujs)2]dw 


(1013) 


-J=S(ujs,  0)el[“st  K^x1  J  exp i[^{-K"(ujs)x  +  ip"(u)s,  0)}(u;  -  ixs)2}du  = 


(1014) 


Ua(t,x)  ~SMJxK„M_rMe 


iuJst—iK (ios )x— i7rsgn  [K" x— 'll)" (ws  ,0)] / 4 


(1015) 


The  phase  is 


<fa(t,  x)  =  1p(u!s,  0)  +  LVst  —  K(us)x  —  7rsgn  [K"x  —  if" (u si  0)]/4  (1016) 


Differentiating  the  phase,  ipa(t  ,x  ),  we  have 


ki(t,x)  =  ~—ipa{t,X ) 
ox 


dip  dus  dus  dK(us)dus1\ 

ft  +  ( +  km 

dus  dus  J  ox 

dip  f  t  +  -0'(^>°)M  duJs  ,  T,,  , 

J7  +  t~x{ - x - j  ~te+K(u,) 


(1017) 

(1018) 


(1019) 


x 


(1020) 


However 


and  hence 


Also, 


which  gives 


A» 


t  +  ip'{u,  0) 
x 


ki(t,x)  =  K(u>s ) 


.  d 

Wi[t,X)  =  —<pa(t,x) 


dt 
dip  dus 


dus  dt 
dip 


dus  dK(us)dus 

+ +c us-x- 


dus 

dip 

duq 


+ 1  —  x 


+ 1  —  X 


dt 
dK  (cns 


du*  dt 


du. 


dus 

dt 


+  W s 


t  +  ip'(u,  0) 


X 


dus 

dt 


+  W s 


Ui(t,x)  =  Us 


(1021) 

(1022) 

(1023) 

(1024) 

(1025) 

(1026) 

(1027) 
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